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ABSTRACT 


An  attempt  is  made  to  find  the  rules  that  correlate 
the  channel  geometry  of  gravel  rivers  and  gravel  canals  with 
discharge  and  bed  material.  Five  gravel  rivers,  which  flow 
out  of  lakes,  were  investigated  a  short  distance  downstream 
of  the  lake  outlets.  Four  laboratory  tests  were  made  to  ex¬ 
tend  the  ranges  of  discharge  and  grain  size.  Use  was  also 
made  of  gravel  river  and  gravel  canal  data  available  in  var¬ 
ious  publications.  The  relation  between  tractive  force  and 
sediment  size,  a  fourth-root  flow  formula  and  a  modified 
version  of  Lacey’s  width-discharge  relation,  were  used  to 
derive  a  set  of  regime  equations  for  gravel  rivers  and  gravel 
canals  with  low  sediment  discharge  and  low  depth.  Methods  of 
sampling  coarse  river  beds  are  discussed  and  it  is  concluded 
that  areal  grid  samples,  worked  out  by  number  in  terms  of  the 
equivalent  spherical  diameter,  or  in  terms  of  the  intermediate 
axis,  can  be  used  as  approximations  to  the  standard  sieve  curve 
of  the  surface  layer  of  a  gravel  paved  channel. 
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CHAPTER  i,  INTRODUCTION 


1 . 1 _THE  PKOPLBM 

Almost  .ill  the  large  rivers  in  Alberta  and  in  British 
Columbia  flow  on  gravel  or  boulder  beds  during  the  major  part 
of  their  course  in  these  two  provinces.  Interferences  such  as 
dams ,  diversions  and  dikes,  are  becoming  more  and  more  frequent 
as  the  development  of  natural  resources  progresses  towards  and 
into  the  mountainous  areas  at  the  headwaters  of  these  rivers. 

To  appreciate  the  consequences  of  interferences,  the  factors 
governing  the  equilibrium  or  regime  of  gravel  rivers  ought  to 
be  known.  It  is  the  purpose  of  this  thesis  to  explore  one 
special  aspect  of  this  problem,  the  equilibrium  for  low  sedi¬ 
ment  transport  rates. 

A  number  of  gravel  rivers  below  natural  lakes  were  studied 
during  the  summer  of  1962.  During  the  following  winter  a  few 
flume  experiments  were  made  to  investigate  the  no-transport 
equilibrium  of  gravel  channels.  Based  on  the  results  of  this 
work  an  attempt  will  be  made  t©  describe  the  characteristic 
properties  of  a  gravel  river  with  low  sediment  discharge.  It 
is  hoped  that  this  study  will  help  in  predicting  the  changes 
that  take  place  in  a  gravel  river  downstream  of  a  newly  con¬ 
structed,  dam. 

1.2  NOTATIONS  AND  DEFINITIONS 

1.2.1  Notation.  All  letter  symbols  used  in  this  thesis  will 
be  defined  where  they  first  appear.  For  convenience,  the  more 
commonly  used  symbols  and  abbreviations  will  be  listed  in 
Appendix  A. 


( 


1.2.2  Sand  , Gr a  ve  1 ,  Rou 1 d  e  r  s .  Material  consisting  mainly  of 

grains  coarser  than  2  nun  will  be  defined  as  gravel.  Sand  will 
mean  material  predominantly  in  the  grain  size  range  from  0.06 
mm  to  2  mm.  It  is  usually  obvious  whether  the  bed  material  of 
a  particular  river  reach  is  either  sand  or  gravel.  River  beds 
consisting  mainly  of  the  intermediate  material  in  the  1.5  to  4 
mm  range  are  uncommon,  (see  Chapter  5). 

No  upper  grain  size  limit  will  be  imposed  on  gravel.  The 
term  ’’boulder"  will  be  reserved  for  single,  exceptionally  large 
stones . 

1.2.3  The  Sediment  Load.  The  part  of  the  total  load  that  moves 
by  sliding,  rolling  or  saltation,  will  be  called  bed  load.  The 
remaining  part  is  carried  in  suspension  and  will  be  called  sus¬ 
pended  load.  The  fraction  of  the  bed  load  that  consists  of  mat¬ 
erial  similar  to  that  of  the  river  bed  will  be  referred  to  as 
bed  material  load,  and  the  remaining  fraction  as  wash  load. 

1.3  _ THE  CONCEPT  OF  REGIME 

1.3.1  Definition.  What  is  meant  here  by  the  synonymous  ex¬ 
pressions  of  "river  in  equilibrium",  or  "river  in  regime",  prob¬ 
ably  needs  some  clarification,  since  there  has  been  considerable 
argument  about  these  concepts,  (Ref.  1,  2,  3,  5). 

Regime  only  applies  to  alluvial  river  reaches.  Such  a 
reach  is  considered  to  be  in  regime  or  in  equilibrium,  when  its 
average  dimensions  such  as  slope,  width,  depth,  meander  wave 
length,  etc.,  do  not  change  appreciably  over  periods  of  many 
years.  The  reach  should  be  a  few  meander  lengths  long.  The 
average  dimensions  may  oscillate  about  a  certain  value,  but  they 
should  not  exhibit  a  definite  trend.  If  meander  cutoffs  occur 
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in  the  reach  under  consideration  the  river  slope  may  oscillate 
about  a  mean  value  with  periods'of  50  or  more  years. 

1.3.2  The  Ideal  Regime  River.  The  idealized,  perfect  regime 
river  has  to  be  visualized  as  flowing  across  a  broad,  very  long 
and  gently  sloping  plane;  built  up  by  the  river  itself  in  earl¬ 
ier  times.  The  downstream  end  of  this  plane  is  formed  by  an 
undest ruct ible  barrier,  across  which  the  river  discharges  at 
supercritical  velocity.  The  state  of  equilibrium  of  any  reach 
along  this  river  will  depend  almost  entirely  on  the  functions 
of  time  that  describe  the  water  and  the  sediment  load  (water 
sediment  complex)  entering  the  reach  at  this  section;  the  pos¬ 
sible  exceptions  to  this  rule  are  factors  such  as  climate  and 
vegetation  along  the  reach. 

Somewhat  simplified  this  can  be  expressed  as: 


spi  ,  SR,  Sec.,  W^,  Lm  =  fns.  (Q,  Wa ,  C,  Sed) 
Where  Sp  Slope  of  the  plane 


1.3.2. 1 


SR  Slope  of  the  river 


average  values 


Sec.  Cross  section  of  the  river 


along  a  reach. 


W 


Width  of  the  meander  belt 


L 


Meander  wave  length 


Q  Discharge 

Wa  Physical  properties  of  water 


C  Sediment  discharge 

Sed  Physical  properties  of  sedi- 


temperature. 


such  as  density,  viscosity, 


Functions  of 


time  . 


ment  such  as  size  distribution, 


resistance  to  abrasion,  etc. 
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The  width  of  the  river  or  the  width  of  the  meander  belt  are  in 
certain  cases  independent  variables;  examples  being  rivers 
meandering  in  restricted  channels  that  were  cut  under  different 
geological  conditions,  flume  experiments,  and  rivers  trained 
between  rigid  banks. 

Equation  1.3. 2.1  is  simply  a  list  of  possible  causes  and 
their  effects.  It  does  not  imply  that  all  the  independent  var¬ 
iables  shown  actually  influence  all  or  any  of  the  dependent 
variables.  Experience  with  unlined  irrigation  canals  has  proven 
conclusively  that  a  relation  similar  to,  but  simpler,  than 
equation  1.3.2. 1  (no  meandering,  independent  variables  essent¬ 
ially  constant  with  time),  is  unique  within  narrow  limits,  (Ref. 
2 , 5 , 6 , 7 , 8 , 9 ) .  Sediment  discharge  was  not  measured  in  most  of 
the  canals  mentioned  in  the  above  references.  The  bed  load  was 
usually  very  small  but  the  suspended  load  of  some  canals  was 
considerable . 

Geologists,  until  recently,  evaded  the  question  of  depend¬ 
ency  by  talking  of  a  graded  river  and  defining  it  as  a  river  in 
which  the  slope  is  adjusted  to  transport,  with  prevailing  chan¬ 
nel  characteristics,  the  sediment  supplied  at  the  head  of  the 
graded  reach.  (Ref.  1.) 

1.3.3  Application  of  the  Regime  Concept  to  Rivers.  Steadily 
changing  geological  and  climatic  conditions  make  the  perfect 
state  of  equilibrium,  as  described  in  1.3.2,  a  rather  unlikely 
occurrence.  Particularly  at  present,  only  a  few  thousand  years 
after  an  ice  age,  ideal  equilibrium  cannot  be  assumed  lightly 
for  the  majority  of  alluvial  river  reaches. 

The  concept  of  the  regime  river  has  become  valuable  to 


engineers  and  geologists  in  spite  of  these  limitations  for  the 
following  reason:  a  river  that  aggrades  or  degrades  usually 
does  so  very  slowly,  or  in  other  words,  the  trends  mentioned 
in  1.3.1  will  be  slow.  On  a  long  time  average  the  sediment  load 
increase  or  decrease,  along  a  reach  of  a  few  meander  lengths, 
is  small  compared  with  the  load  transported  through  that  reach. 
This  probably  accounts  for  the  near  impossibility  of  disting¬ 
uishing  between  aggrading,  degrading  and  regime  rivers  on  the 
basis  of  channel  geometry  and  hydraulic  measurements  not  cover¬ 
ing  periods  of  several  life  times.  (Ref.  4,  p.45-49.)  Summar¬ 
izing,  one  may  say  that  most  rivers  are  not  truly  in  equilibrium 
but  their  channel  geometry  and  hydraulic  characteristics  can  be 
expected  to  resemble  the  equilibrium  state  closely. 

Work  done  by  the  U.S.G.S.  (Ref.  10),  sugge sis  that  equation 
1.3.2. 1  is  applicable  to  sand  rivers,  as  far  as  Q  as  an  inde¬ 
pendent  variable  is  concerned.  For  the  majority  of  rivers 
investigated  the  sediment  discharge  was  not  available. 

For  gravel  rivers  the  uniqueness  of  the  relation  stated  in 
equation  1.3.2. 1  has  never  been  established,  particularly  where 
width  as  a  dependent  variable  is  concerned.  The  equilibrium  may 
depend  mainly  on  the  way  in  which  it  was  reached.  The  geologi¬ 
cal  history  of  such  rivers  is  therefore  important. 

1.3.4  Regime  of  Rivers  below  Lakes.  A  river  without  sediment 
load  will  degrade  its  bed  eventually,  even  if  it  consists  of 
coarse  gravel  that  does  not  move  at  the  highest  floods.  Chemi¬ 
cal  and  physical  agents  will  disintegrate  the  stones.  In  add¬ 
ition,  if  climatic  factors  cause  the  discharge  of  such  a  river 
to  decrease,  it  might  not  be  able  to  change  its  channel  to  suit 
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the  new  conditions.  This  means  that  river  reaches  just  below 
large  lakes  downstream  to  the  first  sediment  supplying  tribut¬ 
ary,  may  not  be  in  regime,  even  if  they  are  alluvial. 
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CHAPTER  2  THE  I-'lELL)  WORK 


Some  general  aspects  of  the  field  work  will  be  described 
in  paragraphs  1  to  8.  In  paragraphs  9  to  15  the  relevant  data 
of  each  measuring  reach  will  be  presented  separately.  The 
important  measurements  of  all  the  reaches  are  summarized  in 
Tables  2  to  9. 

2 . 1  THE  SELECTION  OF  SUITABLE  RIVERS 

For  the  proposed  study  of  the  low- transport  equilibrium, 
a  gravel  river  had  to  meet  a  rather  large  number  of  limiting 
conditions,  which  are  stated  below. 

a)  In  order  to  make  the  assumption  of  a  low  sediment 
transport  equilibrium  possible,  suitable  study  reaches  had  to 
be  on  rivers  that  flow  out  of  lakes,  preferably  not  more  than 
a  few  miles  downstream  of  such  a  lake. 

b)  As  pointed  out  in  1.3.1,  the  study  reaches  had  to  be 
alluvial.  To  have  a  small  tributary  stream  join  the  main  river 
just  upstream  of  the  study  reach  was  desirable,  particularly  if 
the  tributary  had  a  small  discharge  and  a  large  sediment  load. 
(See  1.3.4) 

c)  The  discharge,  at  the  time  when  measurements  were  made, 
had  to  be  known,  and  a  comparison  between  this  flow  and  the  gen¬ 
eral  flow  characteristics  of  the  river  had  to  be  possible.  This 
meant  that  the  rivers  had  to  be  gauged  and  that  flow  records  of 
at  least  twenty  years  had  to  be  available. 

d)  The  ideal  river  would  have  had  a  fairly  constant  dis¬ 
charge  all  year  round,  but  this  could  obviously  not  be  expected. 
Based  on  the  assumption  that  gravel  rivers  with  very  little  sus- 
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pended  sediment  and  small  bed  load  charge  are  morphologically 
inactive  at  low  stages,  it  was  specified  only  that  the  river 
should  have  a  small  ratio  between  floods  of  high,  and  floods  of 
low  frequency.  The  question  of  what  discharge  should  be  con¬ 
sidered  responsible  for  a  given  river  channel  (dominant  dis¬ 
charge),  loses  some  of  its  importance,  if  this  type  of  river  is 
studied . 

e)  The  rivers  had  to  be  reasonably  free  from  human  inter¬ 
ferences  to  ensure  equilibrium  conditions. 

f)  Reasons  of  time  and  money  imposed  some  limitations  on 
the  distances  that  could  be  travelled  to  obtain  field  data. 

g)  Only  river  reaches  accessible  at  least  by  a  four-wheel 
drive  vehicle  could  be  considered.  Financial  reasons,  and  the 
weight  of  the  equipment  excluded  access  by  plane  or  boat. 

With  the  help  of  maps,  and  after  consulting  the  Alberta 
Water  Resources  Branch,  it  was  concluded  that  no  river  in 
Alberta  meets  the  requirements  outlined  above  sufficiently  to 
justify  studying  it. 

The  author  knew  vaguely  of  some  possibilities  in  British 
Columbia.  These  were  confirmed  through  correspondence  with  the 
B.C.  Water  Rights  Branch,  and  with  the  aid  of  aerial  photographs 
A  considerable  number  of  rivers  in  Central  B.C.  south  of  Quesnel 
appeared  to  be  possible.  Finally,  the  Chilko,  Taseko,  Quesnel 
and  Cariboo  Rivers,  below  the  respective  lakes  of  the  same  names 
were  chosen.  Fig.  1  shows  their  general  location.  It  was  left 
to  be  decided  in  the  field,  exactly  what  reaches  would  be  stud¬ 
ied.  Not  all  the  conditions  specified  above  were  satisfied  by 
each  river,  but  they  all  met  the  important  conditions,  2.1.b 
and  d,  alluvial  reach  and  similar  flood  discharges  each  year, 
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All  four  rivers  are  within  one  hundred 
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Lake,  so  that  travel  between  them  did  not 
consume  an  unreasonable  part  of  the  time  spent  in  the  field. 
Storage  dams  have  been  considered  by  various  agencies  on  all 
the  lakes  mentioned  above.  Detailed  maps  of  the  lake  outlet 
areas  were  therefore  available. 

Measurements  made  by  the  B.C.  Water  Rights  Branch  on  the 
Thompson  River  below  Kamloops  Lake  will  also  be  included  in 
this  study.  The  author  spent  only  one  day  there  in  late 
September,  sampling  the  bed  material. 

2 . 2  ORGANIZATION  OF  FIELD  WORK 

Each  river  was  visited  twice;  first  in  June  or  July,  and 
again  in  late  August  or  September.  The  main  objectives  of  the 
first  visit  were  the  selection  of  one  or  two  suitable  measuring 
reaches,  and  the  measurement  of  water  surface  slope,  and  water 
velocity  at  high  stage.  During  the  second  visit  to  a  river  in 
fall  the  bed  material  was  sampled  and  cross  sections  were  meas¬ 
ured  at  most  stations. 

The  field  work  would  have  been  simplified  considerably  if 
the  rivers  could  have  been  visited  first  at  low  stage. 

2 . 3  SELECTION  OF  STUDY  REACHES 

Along  measuring  reaches  the  rivers  had  to  be  reasonably 
straight  and  had  to  have  a  well  defined  channel.  Dense  bush 
and  difficulties  of  access  usually  limited  the  number  of  poss¬ 
ible  measuring  reaches.  For  the  same  reasons  most  of  the 
reaches  that  were  finally  selected  were  not  as  long  as  would 
have  been  desirable  for  obtaining  a  representative  slope.  This 
is  certainly  one  of  the  most  seveie  1 imi tat  ions  of  this  study. 
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Fortunately  three  of  the  six  slope  measurements  could  be  tied 
into  the  much  longer  slope  surveys: made  by  the  International 
Pacific  Salmon  Fisheries  Commission  and  by  the  B.C.  Hydro  and 
Power  Authority.  The  close  agreement  is  reassuring. 

2 . 4  WATER  SURFACE  SLOPE  AND  SURFACE  VELOCITY 

Once  a  reach  was  selected,  three  to  six  stations  were  est¬ 
ablished  along  it  at  fairly  equal  intervals.  These  stations 
were  marked  on  both  sides  of  the  river  with  blazes,  or  with 
flagging,  and  numbered  consecutively  in  a  downstream  direction. 
The  distance  across  the  river  between  corresponding  station  mar¬ 
kers  was  measured  by  triangulation  or  by  stadia.  The  distance 
along  the  river  bank  was  obtained  by  chaining.  On  one  side  of 
the  river^  bench  marks  were  established  at  each  station  and  tied 
into  each  other.  Water  surface  elevations  were  usually  meas¬ 
ured  twice  at  each  station,  once  during  the  first  visit  to  a 
river  and  once  during  the  second  visit  in  fall.  The  water  sur¬ 
face  velocity  was  obtained  at  about  two  locations  along  each 
reach.  It  was  measured  by  timing  a  float  along  a  short  distance. 

2 . 5  SAMPLING  THE  BED  AND  BANK  MATERIAL 

All  the  river  beds  had  a  very  distinct  pavement,  consist¬ 
ing  of  the  coarsest  stones  available  from  the  bed  material. 
Diameters  ranging  from  3  to  15  inches  were  most  commonly  en¬ 
countered.  This  pavement  material  was  sampled  by  picking  stones 
at  random,  and  measuring  their  three  axes.  This  areal  grid 
sampling  technique  was  first  described  by  Wolman  in  Ref.  11,  and 
was  used  extensively  by  Qureshi  (Ref.  12).  Wolman  only  measur¬ 
ed  the  intermediate  axis  of  each  stone.  This  is  justified  if 
all  the  samples  are  of  reasonably  equal  roundness.  However  the 
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samples  that  were  collected  for  this  study  had  widely  varying 
average  shapes.  To  take  this  into  account  all  three  axes  of 
each  stone  were  measured,  as  described  by  Qureshi  in  Ref.  12. 
Random  picking  appears  to  be  the  only  possible  method  of  samp¬ 
ling  very  coarse  material,  if  the  dredging  and  sieving  of 
samples  as  large  as  a  few  cubic  yards  is  excluded  for  financial 
reasons . 

All  the  random  picking  of  stones  from  under  water  was  done 
by  walking  over  the  area  to  be  sampled,  making  long  steps  with 
closed  eyes,  and  after  each  step  picking  up  the  stone  that  hap¬ 
pened  to  be  under  the  toes.  In  this  way  the  sample  area  was 
covered  by  a  zigzag  course.  Usually  the  general  direction  was 
kept  constant.  A  sample  of  60  to  80  stones  covers  then  an  area 
of  about  3  by  100  feet.  Sometimes  a  broader  and  shorter  area 
was  sampled  by  moving  back  and  forth  across  it,  taking  care  not 
to  cover  any  particular  place  more  than  once.  Between  50  and 
80  stones  per  sample  would  have  been  desirable.  Usually  the 
effect  of  the  very  cold  water  resulted  in  smaller  samples. 

Even  at  low  stage  in  fall,  most  of  the  rivers  were  so  swift 
that  it  was  impossible  to  collect  stones  from  a  depth  greater 
than  two  feet.  Fortunately  all  the  rivers,  except  the  Taseko, 
had  very  clear  water.  The  bed  material  could  be  seen  all  across 
the  channel,  either  from  shore  or  from  a  boat,  and  an  estimate 
could  be  made  of  how  well  a  particular  sample  represented  the 
bed  pavement  of  the  central  part  of  a  channel. 

Most  rivers  were  quite  low  during  the  second  visit  in 
September.  Often  places  could  be  found  above  water,  where  the 
pavement  appeared  to  be  very  similar  to  that  existing  in  the 
main  channel.  Tpis  was  always  a  welcome  sight  as  the  discomfort 
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of  picking  stones  from  two  feet  of  cold,  swift,  water  could 
then  be  avoided.  Random  sampling  above  water  was  done  by 
stretching  a  survey  tape  across  the  sampling  area  and  picking 
up  stones  below  every  foot-marker.  A  f  ew  samples  were  col  lected 
by  picking  and  measuring  all  the  stones  larger  than  some  mini¬ 
mum  value,  which  were  exposed  within  a  certain  small  area. 

The  volume  and  equivalent  spherical  diameter  (diameter  of 
sphere  of  equal  weight),  were  calculated,  using  the  formulas 
given  in  Ref.  12. 

V  =  cz  (  a  b  )z  2.5.1 


<t>  =  ({cl  a^  bJ  )3  2.5.2 

Where  a,  b,  and  c,  are  the  largest,  intermediate,  and  smallest 
axis  of  the  stone;  V  is  the  volume  and  <^>  is  the  diameter  of  a 
sphere  of  equal  weight.  Fig.  2  shows  a  useful  nomogram  for  cal¬ 
culating  V  and  4>  .  The  rather  unusual  selection  of  units  re¬ 
sults  from  the  fact  that  a,  b,  and  c,  were  measured  with  a 
piece  of  survey  tape,  but  it  was  more  convenient  to  have  in 
inches,  to  facilitate  comparison  with  other  studies. 

Equation  2.5.1  could  only  be  verified  in  a  few  cases, 
where  the  stones  of  a  sample  were  small  enough  to  be  weighed  on 
a  small  balance.  The  results  are  shown  in  Table  1  together  with 
some  data  from  Ref.  14.  The  weight  was  calculated  for  a  speci¬ 
fic  gravity  of  2.65.  Agreement  appears  reasonably  good.  Con¬ 
sidering  that  <£>  only  depends  on  the  third  root  of  V,  one  can 
expect  4*  to  be  accurate  to  within  10  percent.  It  is  probably 
a  few  percent  too  small,  but  the  scanty  data  did  not  justify 
the  use  of  a  correction  coefficient  in  equation  2.3.1. 

Frequency  curves  (Log  vs.  ’’percent  by  number  finer  than”) 


were  plotted  for  each  sample.  These  curves  are  included  in  this 
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report  (see  2.9  to  2.15).  The  more  important  gravel  samples 
were  also  plotted  as  "logarithm  of  intermediate  axis  vs.  per¬ 
cent  by  number  finer  than",  as  suggested  by  Wolman.  Only  one 
of  these  curves  is  shown  here  (Fig.  49),  but  the  results  of 
analysing  the  samples  in  this  fashion  can  be  seen  in  Table  2. 

The  sampling  technique  described  above  is  comparatively 
new  and  not  generally  accepted.  It  will  therefore  be  discussed 
briefly  in  Chapter  4. 

At  a  few  favourable  places,  such  as  road  cuts,  samples  of 
the  bank  material  were  taken  with  a  shovel.  This  material  was 
split  and  sieved  wet,  down  to  mesh  number  four.  The  finer 
fraction  was  taken  to  Edmonton,  and  sieved  to  mesh  200.  The 
methods  used  for  sampling  and  wet  sieving  were  rather  crude,  so 
that  the  resulting  sieve  curves  may  not  be  representatives-  "  v' . 


2.6  CROSS  SECTIONS 


At  almost  all  s 
with  level  and  rod, 
possible.  The  cross 
above  bankfull  stage 
sections  are  complet 
could  be  used  for  ec 
spans  the  river.  Th 
estimated,  using  fie 
shape  of  the  main  ri 
and,  where  feasible, 
applied  to  the  water 


tations cross  section  measurements  were  made 
and  by  wading  out  into  the  river  as  far  as 
sections  usually  extend  from  a  point  well 
toasimilar  point  across  the  river.  The 
e  only  at  those  stations  where  the  boat 
ho  sounding,  or  where  a  bridge  or  cable 
e  missing  parts  of  the  cross  sections  were 
Id  book  entries  on  the  apparent  depth  and 
ver  channels,  measured  surface  velocities, 
the  Manning  equation.  This  equation, 
surface  slope  measurements  at  high  and  low 


stage,  gives, under  suitable  circumstances,  a  good  estimate  of  the 
missing  water  depth  data. 
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The  sketch  below  shows  a  typical  river  cross  section.  The 
dashed  line  indicates  the  missing  parts. 


Index  1  refers  to 
high  stage,  index 
2  to  low  stage. 


Manning ' s 

Q,  = 

Qz  = 
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and  rearranged: 
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2.6.1 


All  the  terms  in  equation  2.6.1,  expect  d,  ,  are  known.  As  no 
width  to  depth  ratios  of  less  than  15  were  encountered,  hydrau¬ 
lic  radius  and  wetted  perimeter  were  usually  replaced  by  an 
average  width,  and  by  the  depth  of  the  central  part  of  the 
channel,  so  that  the  product  width  times  depth  is  equal  to  the 
section  area.  Depth  estimates  with  equation  2.6.1  were  made 
for  known  sections  to  test  this  method.  Very  close  agreement 

j 

was  reached,  whenah  was  larger  than  about  and  when  the  sec¬ 
tion  was  trapezoidal,  and  in  a  straight  reach.  (See  2.13.) 

It  should  be  kept  in  mind  that  some  of  the  estimated  sections 
are  based  on  scanty  data,  particularly  the  Quesnel  and  Taseko 
sect  ions . 

In  future  d  will  stand  for  the  average  channel  depth  in 
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the  central  part  of  the  river,  as  shown  above.  The  gravel 
rivers  studied  were  found  to  have  fairly  flat  beds  over  a  major 
part  of  the  width,  so  that  d  is  usually  well  defined, 
d*  will  stand  for 


where  Ws  is  the  water  surface  width  and  A  the  cross  section 
area.  In  shallow  wide  rivers,  d*  is  only  a  few  percent  larger 
than  the  hydraulic  radius. 

2 . 7  DISCHARGE  RECORDS 

With  the  exception  of  the  Taseko,  all  the  rivers  that  were 
studied  had  gauging  stations  very  near  the  measuring  reach. 

Flood  frequency  curves  (annual  maximum  daily  mean  discharge  vs. 
frequency  of  occurrence)  were  plotted  on  logarithmic  probability 
paper.  (Fig.  3).  The  same  data  were  also  plotted  on  arithmetic 
probability  paper.  Tfte  logarithmic  plots  are  given  here,  be¬ 
cause  their  fitting  lines  are  considerably  straighter. 


2.8  REGIME  CONSIDERATIONS 

In  order  to  facilitate  the  comparison  of  the  equilibrium 
state  in  different  rivers,  a  certain  typical  discharge  has  to 
be  assigned  to  each  river.  What  discharge  can  be  considered  as 
typical  poses  a  difficult  question.  Possibly  there  is  no  com¬ 
pletely  correct  answer,  as  the  varying  discharges  of  a  river 
might  produce  a  channel  that  could  not  have  been  obtained  by  a 
certain  constant  flow.  In  Ref.  10  the  problem  was  solved  by 
using  discharges  of  equal  frequency  of  occurrence  on  all  the 
rivers,  without  trying  to  decide  whether  these  discharges  actu¬ 
ally  could  have  formed  the  channel.  In  this  study  more  is 


known  about  the  recent  geological  history  of  the  measuring 
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reaches  than  in  Ref.  10,  therefore  an  attempt  has  been  made  to 
select  representative  or  dominant  discharges  that  could  have 
formed  the  particular  channel. 

Some  of  the  study  reaches  are  very  close  to  lake  outlets. 
Almost  no  sediment  is  supplied  to  the  river  upstream  of  these 
rea'ches,  so  that  the  bed  and  bank  erosion  and  the  flattening 
of  the  slope  caused  by  extreme  floods  cannot  be  repaired.  The 
typical  or  dominant  discharge  is  therefore  a  flood  of  low  fre¬ 
quency.  As  pointed  out  in  1.3.4,  there  is  no  guarantee  that 
equilibrium  exists  at  all,  even  during  an  extremely  high  flood. 

Other  measuring  reaches  are  located  several  miles  down¬ 
stream  of  lakes.  Tributaries  supply  sediment  to  the  river 
between  the  lake  outlet  and  the  study  reaches.  Here  the  eros¬ 
ion  caused  by  an  extreme  f lood  can  to  some  extent  be  repaired 
by  the  river.  An  extreme  flood  lasting  perhaps  a  few  weeks 
would  also  not  be  able  to  flatten  the  river  slope.  At  dominant 
discharge  such  reaches  must  just  be  able  to  transport  the 
coarsest  fractions  of  sediment  entering  from  tributaries. 

Floods  with  a  recurrence  interval  of  from  2  to  5  years  were 
considered  to  be  representative  of  these  reaches.  Whenever 
feasible^ the  highest  discharge  at  which  measurements  had  been 
made  was  taken  as  the  dominant  discharge. 

Instead  of  using  a  dominant  discharge  estimated  on  the 
basis  of  frequency  considerations,  some  people  prefer  to  work 
in  terms  of  a  ”bankfull  discharge”,  determined  from  cross  sec¬ 
tions  and  gauge-discharge  relations.  A  number  of  good  reasons 
can  be  put  forward  for  doing  so,  particularly  if  the  river  is 
meandering  in  a  valley  with  distinct  flood  plains.  The  author 
does  not  believe  that  the  concept  of  ”bankfull  discharge”  is 
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useful  when  dealing  with  incised  stable  channels  shortly  below 
lakes  (e.g.  Thompson  River  below  Kamloops  Lake,  Cariboo  River 
below  Cariboo  Lake).  Presumably  bankfull  discharge  for  these 
channels  would  have  to  be  taken  as  the  discharge  at  which  the 
water  reaches  the  apparent  high  water  line,  as  shown  by  soil, 
vegetation  and  driftwood.  One  wonders  whether  the  elevation  of 
this  line  is  not  primarily  a  function  of  the  climate,  the  type 
of  vegetation  along  the  river  and  the  discharge  regime  of  the 
river.  However,  the  author's  estimate  of  dominant  discharge  is 
almost  identical  with  bankfull  discharge  for  reaches  that  lie 
several  miles  below  lakes  such  as  "Chilko  River  at  Henry's 
Crossing",  or  "Quesnel  River  at  Lawless  Creek",  etc.  Bankfull 
discharge  will  be  given  separately,  wherever  it  differs  apprec¬ 
iably  from  the  assumed  dominant  discharge. 

An  estimate  of  the  typical  cross  section  of  each  measuring 
reach  at  dominant  discharge  was  also  made.  These  estimates  are 
based  mainly  on  the  cross  section  measurements.  On  occasions, 
an  inspection  of  the  river  both  upstream  and  downstream  of  the 
measuring  reach  has  also  been  taken  into  consideration. 

2 . 9  THE  QUESNEL  RIVER  AT  LAWLESS  CREEK 

Description  of  the  River  Upstream  of  the  Measuring  Reach . 

The  Quesnel  River  originates  in  Quesnel  Lake  (area:  105  sq.  m.) 
on  the  western  slope  of  the  Cariboo  Mountains.  The  Dominion 
Water  Resources  gauging  station,  8KH1,  (Quesnel  River  at  Likely, 
B.C.),  is  situated  just  below  the  lake.  At  low  stages  the 
remnants  of  a  small  dam  can  be  seen  in  the  river  near  the  gaug¬ 
ing  station.  This  dam  was  built  in  1898  and  used  for  drying  up 
the  river  in  winter  so  that  gold  could  be  mined  from  its  channel. 
The  high  stages  of  the  Quesnel  River  were  never  greatly  influ- 
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enced  by  the  dam.  About  half  a  mile  below  the  gauging  station 
the  river  enters  the  Red  Canyon.  The  reach  between  the  remains 
of  the  dam  and  Red  Canyon  is  not  suitable  for  measurements  as 
the  river  bed  there  consists  of  angular  boulders  up  to  several 
feet  in  diameter. 

The  Red  Canyon  is  about  5  miles  long.  Huge  boulders  or 
possibly  bedrock  are  exposed  in  the  river.  The  tailings  of  a 
placer  mine,  (Bullion  Mine),  were  once  discharged  into  the  river 
at  a  point  approximately  half  way  down  Red  Canyon.  The  Bullion 
Mine  was  started  in  the  19th  century.  Until  1922  it  operated 
intermittently,  but  from  1922  to  1942  it  was  in  full  operation, 
discharging  around  800,000  cu.  yards  of  sand  and  gravel  into  the 
river  each  year.  This  amounts  to  a  sediment  charge  of  about 
5-20  parts  per  million  during  the  summer  months.  Since  1942  the 
mine  has  been  closed  down. 

Flow  Records.  The  general  flow  characteristics  of  the 
Quesnel  River  can  be  seen  from  Fig.  3  (flood  frequency),  Fig,  6 
(daily  discharges  of  a  few  years),  and  Fig.  7  (flow  duration). 

All  these  plots  are  based  on  the  38  years  of  records  for  the 
gauging  station  8KH1  at  Likely,  B.  C. 

The  Measuring  Reach.  A  measuring  reach  was  selected  at  the 
downstream  end  of  Red  Canyon.  There  the  river  flows  through  a 
circular,  flat  area  which  is  almost  completely  surrounded  by  steep 
cliffs,  (Fig.  4  and  5).  Pictures  1  and  2  were  taken  from  a 
position  high  on  these  cliffs. 

The  river  has  recently  been  degrading  along  the  measuring 
reach.  Comparison  between  Picture  1  and  the  air  photos  taken 
in  1955  (Fig.  5),  illustrates  this.  The  air  photos  were  taken 
at  a  discharge  of  10,000  cfs.  and  show  practically  the  whole 
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gravel  flats  flooded;  yet  on  picture  1  a  discharge  of  14,800  cf,s 
produces  no  spill.  The  river  has  gone  down  about  4  feet  in  7 
years.  This  degradation  appears  to  be  caused  by  the  suspension 
of  operations  at  the  Bullion  Mine  in  1942.  The  behaviour  of  the 
Dominion  Water  Resources  gauge  8KH3  on  the  Cariboo  River  near 
Quesnel  Forks  supports  this  conclusionj  (See  2.10.). 

Cross  Sections  and  Slope.  The  cross  sectional  data  are 
very  scanty^  (Fig.  8),  due  to  two  unforeseen  circumstances.  The 
river  could  not  be  crossed  by  boat,  even  in  fall,  and  there  was 
no  other  feasible  access  to  the  left  side  of  the  measuring  reach. 
The  water  surface  width  was  measured  by  triangulation,  so  that 
the  estimated  position  of  the  left  banks  and  of  the  bar  in  sec¬ 
tion  4.  should  be  accurate  to  within  a  few  feet.  The  sections 
appeared  trapezoidal.  Depth  at  each  section  was  estimated  in 
the  field  and  these  estimates  were  checked  with  equation  2.6.1. 
Fig.  9  shows  the  water  surface  slopes  as  measured  in  July  and  in 
September.  Approximate  energy  lines  are  also  plotted.  A  sum¬ 
mary  of  the  channel  geometry  is  given  in  Table  3. 

Bed  and  Bank  Material.  The  results  of  the  bed  and  bank 
material  sampling  are  listed  in  Table  2a  and  are  plotted  on 
Fig.  10.  Picture  3  shows  the  area  of  sample  4.  The  bed  mater¬ 
ial  in  the  central  part  of  the  channel  is  fairly  well  represented 
by  sample  1.  Samples  4  and  5  were  taken  in  an  old  channel  that 
no  longer  carries  water,  even  at  high  river  stage.  Below  sta¬ 
tion  4  the  bed  material  appears  to  become  somewhat  finer.  All 
the  samples  contain  predominantly  very  well  rounded  material. 

Regime  Considerations.  It  may  seem  paradoxical  in  a  study 
of  the  equilibrium  of  gravel  rivers,  to  use  measurements  taken 
on  a  river  that  is  known  to  be  degrading.  In  the  author's 
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opinion  the  following  points  justify  the  assumption  that  this 
measuring  reach  closely  resembles  a  state  of  equilibrium: 

(a)  The  river  has  formed  its  channel  by  digging  into  its  own 
deposits . 

(b)  An  estimate  of  the  volumes  involved  will  show  that  at  high 
stage  the  bed  load  of  actual  bed  material,  (bed  material  load), 
must  be  smaller  than  1  ppm.  At  low  stage  no  bed  material  tran¬ 
sport  takes  place. 

(c)  The  laboratory  work, (Chapter  3),  showed  how  rapidly  a 
state  resembling  equilibrium  becomes  established  in  the  upper 
part  of  a  flume  after  the  feeding  of  sediment  is  suddenly 
stopped . 

(d)  The  argument  put  forward  under  1.3.3  applies  to  this  reach. 

The  dominant  discharge  will  be  taken  as  14,800  cfs.  Slope 
and  section  measurements  have  been  made  at  this  discharge.  No 
bed  load  transport  could  be  heard  while  doing  it.  This  dis¬ 
charge  corresponds  to  a  3  year  flood. 

The  estimate  of  an  average  section  is  given  on  Table  3. 

It  is  based  on  sections  1  to  4,  because  some  distance  below  the 
measuring  reach  the  river  flows  again  in  a  channel  with  a  simi¬ 
lar  cross  section.  It  is  assumed  that  the  wide  channel  at 
sections  5  and  6  is  caused  by  the  much  finer  bank  material  ? 
(sample  3).  Bank  erosion  in  this  part  of  the  reach  is  shown  on 
Fig.  8,  section  6.  This  is  the  only  erosion  that  could  be  de¬ 
tected  along  the  measuring  reach. 
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2.10  THE  CARIBOO  RIVER  NEAR  QUESNEL  FORKS 

Description  of  the  River.  Until  recently  the  Cariboo  River 
was  called  Quesnel  North  Fork.  It  originates  in  Isaac  Lake  (13 
sq.  m.),  on  the  western  slope  of  the  Cariboo  Mountains.  After 
a  course  of  approximately  60  miles  it  joins  the  Quesnel  River 
at  Quesnel  Forks.  On  the  way  it  flows  through  Lanezi  Lake  (6 
sq .  m.),  Sandy  Lake  (4  sq .  m.),  and  Cariboo  Lake  (4  sq .  m.). 

From  approximately  2  miles  below  Cariboo  Lake,  to  a  point  1  mile 
upstream  of  Quesnel  Forks,  the  river  is  deeply  incised. 

Measurements  were  made  along  two  reaches.  The  upper  one  at 
the  outlet  of  Cariboo  Lake  will  be  described  in  2.11. 

The  Measuring  Reach.  Fig.  4  and  5  show  the  lower  study 
reach  near  Quesnel  Forks.  It  was  selected  because  the  river 
there  meets  the  requirements  of  2.1  reasonably  well.  The  major 
limitations  of  this  reach  are: 

-  It  is  very  short. 

-  It  is  not  straight. 

-  Since  it  is  20  miles  below  Cariboo  Lake,  sediment 
transport  at  high  stage  is  considerably  greater  than 
in  the  other  reaches. 

-  At  low  stage  the  water  surface  slope  is  not  uniform. 

Near  station  3  the  main  flow  shifts  from  the  right 
side  to  the  left  side.  The  channel  is  wider  and 
shallower  there  and  rapids  form  at  low  stage. 

The  pictures  4  and  5  were  taken  along  this  reach. 

Flow  Records.  The  Dominion  Water  Resources  gauging  station 
8KH3  (Cariboo  River  near  Hydraulic)  is  right  on  the  measuring 
reach.  The  gauge  is  located  at  station  5,  and  the  measuring  sec¬ 
tion  is  identical  with  cross  section  1.  Fig.  3  (flood  frequency) 
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and  Fig.  13  (flow  duration),  have  been  prepared  from  the  36 
years  of  records  available  for  this  station.  Fig.  11  and  Fig. 

12  show  some  daily  mean  discharges. 

Sections  and  Slope.  The  cross  sections  are  plotted  on  Fig. 
14.  At  station  1  the  cross  section  below  water  is  based  on  re¬ 
cent  soundings  by  the  Dominion  Water  Resources  Branch.  At  station 
4  the  cross  section  has  been  obtained  by  sounding  off  an  old  cable 
car.  The  slope  measurements  are  shown  on  Fig.  15.  The  river 
slope  at  high  stage  is  probably  best  represented  by  the  energy 
line  st .  1  to  st .  4  of  June  29,  1962.  A  summary  of  the  Hydraulic 
Measurements  is  given  in  Table  4. 

Bed  and  Bank  Material.  The  gravel  samples  are  summarized 
in  Table  2b  and  are  plotted  in  Fig.  16.  The  pictures  6  and  7 
give  some  idea  of  the  appearance  of  the  river  bed  and  banks. 
Samples  2  and  5  represent  the  bed  pavement.  It  appeared  that 
the  majority  of  the  stones  that  were  measured  for  these  two  sam¬ 
ples  were  moved  very  recently.  There  was  no  growth  of  algae  on 
them.  Only  a  few  stones  were  as  well  rounded  as  the  material 
that  had  been  sampled  along  the  Quesnel  River.  Samples  1,  3  and 
4b  represent  the  pavement  along  a  line  3  to  4  feet  below  bank- 
full  stage.  Sample  4a  was  taken  about  1  foot  below  bankfull 
stage.  Fig.  17  shows  the  sieve  curves  of  two  samples  taken 
along  this  reach. 

Regime  Considerations.  An  average  cross  section  of  the 
Cariboo  River  in  this  area  is  shown  ob  Table  4.  It  is  derived 
f  rom : 

(a)  Cross  section  measurements. 

(b)  An  inspection  of  the  river  above  the  measuring  reach. 

(c)  Aerial  photographs. 
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The  dominant  discharge  will  be  taken  as  12,000  cfs  based 
on  the  following  reasons: 

(a)  Slope  and  sections  have  been  measured  at  12,000  cfs. 

(b)  The  two  year  flood  is  only  700  cfs  larger. 

(c)  Some  bed  material  was  probably  in  motion  at  12,000 
cfs.  The  sound  of  rolling  stones  could  be  heard  at 
intervals  of  a  few  seconds. 

Specif ic  Gauge .  Specific  Gauge  curves  were  plotted  for  the 
Water  Resources  gauge  8KH3  at  the  downstream  end  of  the  measur¬ 
ing  reach.  From  1926,  when  the  gauge  was  installed,  until  1935, 
no  change  took  place.  From  1935  to  1949  specific  gauges  went 
up  steadily.  The  1949  stage-discharge  curve  is  about  1  foot 
above  the  original  one  over  a  range  of  discharge  from  5000  cfs 
to  12,000  cfs.  After  1949  specific  gauges  went  down  again  and 
appear  to  have  stabilized  at  about  0.2  feet  below  the  1949 
stage-discharge  curve. 

2.11  THE  CARIBOO  RIVER  AT  THE  OUTLET  OF  CARIBOO  LAKE 

A  general  description  of  the  Cariboo  River  is  given  in 

2.10. 

The  Measuring  Reach.  Fig.  18  and  Fig.  19  show  the  upper 
study  reach  on  the  Cariboo  River.  Rollie  Creek  has  pushed  the 
river  to  the  left  side  of  the  valley.  At  present  Rollie  Creek 
flows  into  Cariboo  Lake.  The  area  marked  Rollie  Creek  Delta  on 
Fig.  19  is  densely  wooded.  No  recent  channel  of  the  creek  could 
be  found  there,  but  the  air  photos  and  old  gold  diggings  show 
that  the  area  belongs  to  its  delta.  Pictures  8  and  9  were 
taken  along  this  reach. 

Flow  Records.  The  discharge  measuring  section  of  the 


Dominion  Water  Resources  gauging  station  8KH13  lies  within  the 
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study  reach.  The  gauge  is  situated  on  Cariboo  Lake  near 
Keithley  Creek.  It  was  installed  in  1961.  Fig.  12  shows  the 
1962  record  of  this  station  together  with  the  record  of  the 
gauging  station  8KII3  near  Quesnel  Forks.  From  a  comparison  of 
all  the  peak  flows  that  occurred  during  the  summer  of  1962,  it 
appears  that  those  measured  at  the  outlet  of  Cariboo  Lake  are 
about  12%  smaller  than  the  corresponding  ones  at  Quesnel  Forks. 
With  this  correction  factor  the  flood  frequency  curve  of  Fig.  3 
can  be  used  for  this  reach. 

The  highest  discharges  of  Rollie  Creek  are  estimated  to  be 
in  the  order  of  500  cfs. 

Sections  and  Slope.  The  cross  sections  are  shown  on  Fig. 
20.  Section  2  is  the  Water  Resources  discharge  measuring  sec¬ 
tion.  The  actual  river  only  begins  approximately  200  feet  fur¬ 
ther  downstream.  Between  section  2  and  section  3  the  river  is 
shallow  for  a  short  distance.  At  low  stage  some  white  water 
forms.  No  bedrock  could  be  discovered  in  the  river,  but  the 
deepest  parts  of  the  channel  could  never  be  seen.  Energy  line 
and  water  surface  slope  are  given  on  Fig.  21.  The  flow  veloc¬ 
ity  at  station  1  is  less  than  0.5  ft  s  .  The  water  elevations 
at  this  station  can,  for  all  practical  purposes,  be  considered 
lake  elevations.  The  water  surface  slope,  st.  2  to  st .  7  of 
July  10,  is  probably  most  representative  of  the  river  slope  at 
high  stage.  A  summary  of  the  section  and  slope  measurements 
is  given  in  Table  5. 

Bed  and  Bank  Material.  The  gravel  samples  are  summarized 
in  Table  2c  and  are  plotted  on  Fig.  23.  The  pictures  10  and  11 
will  give  some  idea  of  the  shape  and  arrangement  of  the  bed 
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pavement.  Sample  5  represents  the  central  parts  of  the  bed; 
samples  3  and  4  the  gently  sloping  left  side  of  the  channel; 
and  samples  1  and  2  the  right  side,  about  3  feet  below  bank- 
full  stage.  The  stones  of  all  the  samples  were  covered  with 
algae  and  seemed  to  have  been  at  rest  for  many  years.  At  first 
glance  the  stones  on  the  right  side  appeared  well  rounded,  but 
on  closer  inspection  approximately  half  were  found  to  be  split. 
This  is  probably  a  result  of  freeze-thaw  action,  (picture  11). 

On  the  left  side  rounded,  flat  and  angular  stones  were  mixed. 

The  author  has  no  explanation  for  the  presence  of  the  angular 
stones . 

Fig.  22  shows  three  sieve  curves  of  material  taken  from 
the  Rollie  Creek  delta. 

Some  Regime  Considerations.  The  argument  put  forward  in 
paragraph  1.3.4  applies  to  this  reach.  It  lies  immediately 
below  a  lake  and  no  sediment  load  is  supplied  to  it  at  the  lake 
outlet.  The  river  is  therefore  unable  to  repair  any  erosion 
caused  by  extreme  floods.  A  state  of  equilibrium  could  only 
exist  during  such  a  high  flood.  The  dominant  discharge  is  taken 
as  17,000  cfs,  which  corresponds  approximately  to  the  fifty  year 
flood.  The  maximum  discharge  on  record  is  20,000  cfs  in  Quesnel 
Forks,  which  corresponds  to  around  17,000  cfs  at  this  measuring 
reach.  The  dimensions  of  an  average  section  at  dominant  dis¬ 
charge  are  given  in  Table  5.  Depth  has  been  obtained  by  extra¬ 
polating  the  stage  discharge  curve  of  the  Water  Resources  gaug¬ 
ing  station. 

2.12  THE  TASEKO  RIVER  BELOW  TASEKO  LAKE 

Description  of  the  River.  The  Taseko  River  has  its  origin 
on  the  eastern  slope  of  the  Coast  Range  of  British  Columbia. 
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After  flowing  in  a  NW  direction  for  about  10  miles  it  enters 
Taseko  Lake  (total  area  13  sq .  in.).  The  delta  of  the  Tchaikazan 
River  splits  Taseko  Lake  into  two  parts,  the  upper  and  the  lower 
lake  . 

After  leaving  the  lower  Taseko  Lake,  the  river  takes  a 
northerly  course  across  the  semi-arid  Chilcotin  Plateau.  Fifty- 
five  river  miles  downstream  of  the  lake  outlet  the  Taseko  River 
joins  the  Chilko  River.  In  summer  the  discharge  of  the  Taseko 
River  is  almost  all  glacial  melt  water,  carrying  a  very  fine, 
whitish  silt  load.  The  Taseko  Lakes  are  not  deep  enough  for  this 
silt  to  settle  out.  Therefore,  even  downstream  of  the  lake,  one 
cannot  recognize  any  features  of  the  river  bed  through  more  than 
approximately  1  foot  of  water.  Taseko,  incidentally,  means 
"white  water"  in  the  language  of  the  Chilcotin  Indians. 

The  Measuring  Reach.  A  short  reach  was  studied  on  the  Tas¬ 
eko  River  one  mile  downstream  of  the  outlet  of  Taseko  Lake.  Fig. 
24  is  a  sketch  of  the  area  and  on  Fig.  25  the  lake  outlet  can  be 
seen  in  stereo.  The  situation  is  similar  to  that  at  the  outlet 
of  Cariboo  Lake.  The  delta  of  Beece  Creek  has  pushed  the  Taseko 
River  to  the  left  side  of  the  valley.  On  this  delta  Beece  Creek 
splits  up  into  many  channels.  Some  discharge  into  the  river, 
others  into  the  lake.  In  the  summer  of  1962  almost  all  the 
channels  carried  water,  and  some  new  ones  were  in  the  process 
of  being  formed.  The  larger  channels  have  a  bed  of  coarse 
gravel  with  some  boulders. 

At  the  upper  end  of  the  measuring  reach  the  river  is  some¬ 
what  obstructed  by  large  boulders  that  either  fell  off  the 
mountainside  to  the  left,  or  were  brought  down  by  avalanches. 

An  old  cable  car  crossed  the  river  in  the  centre  of  the  measur¬ 
ing  reach,  but  the  installation  collapsed  when  an  attempt  was 
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made  to  use  it  for  an  echo  sounder  traverse. 

The  measuring  reach  is  not  typical  of  the  Taseko  River 
in  this  area,  due  to  the  presence  of  a  long  bar  on  the  left 
side.  However  dense  bush  made  it  necessary  to  take  advantage 
of  this  rather  unusual  bar  for  levelling  and  chaining. 

Pictures  12-14  show  this  reach. 

Flow  Records.  Only  intermittent  flow  records  of  the 
Taseko  River  are  available.  They  are  for  the  summers  of  1929 
and  1930,  and  were  measured  at  Nemaia  Crossing,  15  miles  down¬ 
stream  of  the  lake  outlet.  There  are  also  about  30  discharge 
measurements  available.  They  were  made  between  1929  and  1949, 
both  at  Nemaia  Crossing  and  at  the  lake  outlet.  An  estimate  of 
the  Taseko  discharge  can  be  made  in  the  following  manner.  The 
Chilko  River  is  gauged  at  the  outlet  of  Chilko  Lake  (8MA2),  and 
forty  miles  further  downstream  just  below  the  confluence  with 
the  Taseko  River  (8MA1).  After  about  mid-June,  the  summer-dry 
Chilcotin  Plateau  contributes  only  slightly  to  the  discharge  of 
the  Taseko  and  the  Chilko  River.  For  the  summer  months,  one 
can  therefore  estimate  the  discharge  of  the  Taseko  River  at 
Taseko  Lake  from  the  difference  between  the  two  Chilko  records. 
Comparison  with  the  existing  Taseko  records  indicates  that  by 
taking  93%  of  the  difference  of  the  two  Chilko  records,  a  good 
estimate  is  obtained  of  the  Taseko  River  discharge  at  the  lake 
outlet.  No  estimate  of  the  Taseko  River  discharge  can  be  made 
for  May  or  for  the  first  half  of  June.  This  is  no  great  handi¬ 
cap,  since  the  existing  Taseko  records  and  information  obtained 
locally^ indicate  that  Taseko  Lake  usually  reaches  its  highest 
stage  only  in  July  or  August.  Knowing  that  the  lake  is  4,300 
feet  above  sea  level  and  that  the  drainage  basin  is  heavily 
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glaciated,  this  is  not  surprising. 

The  flood  frequency  curve  of  Fig.  3  and  the  daily  discharges 
plotted  on  Fig.  26,  are  based  on  the  two  Chilko  records.  The 
correction  factor  mentioned  above  has  not  been  applied  to  them. 

The  discharge  of  Beece  Creek  has  been  measured  occasionally. 
A  maximum  of  520  cfs  is  on  record. 

Section  and  Slope.  Cross  sections  are  plotted  on  Fig.  27. 
Depth  has  been  estimated  with  equations  2.6.1  and  from  field 
notes.  It  was  possible  by  wading  out  from  the  right  bank  to 
obtain  depth  measurements  in  the  main  current.  It  is  unlikely 
that  the  sections  are  much  deeper  further  out.  Air  photos  and 
inspection  of  the  river  indicate  that  only  section  2  can  be  con¬ 
sidered  as  typical  of  the  Taseko  River  just  below  the  lake.  Its 
dimensions  for  August  21,  1962  are  listed  in  Table  6. 

The  slope  measurements  are  shown  on  Fig.  28,  together  with 
the  river  slope  as  it  appears  on  a  large  scale  map  of  the  B.C. 
Hydro  and  Power  Authority.  The  rapids  near  the  large  boulders 
(picture  14),  appear  as  the  steep  drop  at  the  upstream  end  of 
the  slope  plot.  Station  0c  is  well  above  these  rapids.  The 
water  surface  at  station  3  is  high  because  the  water  piles  up 
against  the  upstream  end  of  the  bar.  The  water  surface  slopes 
st.  1  to  st.  5  of  August  13  and  21  are  probably  representative 
of  the  river’s  slope. 

Bed  and  Bank  Material.  A  summary  of  the  bed  material 
samples  is  given  in  Table  2d  and  the  same  samples  are  plotted 
on  Fig.  29.  The  pictures  12  -  14  show  some  of  the  areas 
that  were  sampled.  Samples  2,  5  and  6,  are  believed  to 
represent  the  bed  material  of  the  deeper  parts  of  the  channel. 
Sample  3  represents  the  shallow  left  side  of  the  channel  and 
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sample  1,  the  steep  banks  of  the  right  side.  The  majority  of 
the  stones  were  well  rounded. 

Two  samples  were  collected  with  a  shovel  on  the  right  bank. 
The  sieve  curves  are  shown  on  Fig.  30.  These  samples  should 
give  some  idea  of  the  material  that  Beece  Creek  brings  into  the 
Taseko  River. 

Regime  Considerations.  This  reach  lies  only  a  short  dis¬ 
tance  below  a  lake.  The  sediment  load,  particularly  the  bed 
material  load,  must  be  extremely  small.  Since  bed  erosion  caused 
by  extreme  floods  can  only  be  repaired  slowly,  the  assumption  of 
a  fairly  high  dominant  discharge  seems  justified.  It  will  be 
taken  as  6,400  cfs,  which  corresponds  to  the  once  in  ten  years 
flood.  A  discharge  of  5,400  cfs  has  been  measured  at  Nemaia 
Cr  ossi ng . 

The  dimensions  of  an  average  section  are  given  in  Table  6. 


2.13  THE  CHILKO  RIVER  AT  HENRY’ S  CROSSING 

Description  of  the  River.  The  Chilko  River  originates  in 
Chilko  Lake  (area  70  sq.  m.)  on  the  eastern  slope  of  the  Coast 
mountains  of  British  Columbia.  From  Chilko  Lake  the  river  flows 
in  a  north  westerly  direction  across  the  Chilcotin  plateau.  After 
a  course  of  approximately  50  miles,  it  joins  the  Chilcotin  River 
near  Redstone,  B.C.  Fig.  31  is  a  sketch  of  the  first  nine  miles 
of  the  Chilko  River,  showing  the  two  reaches  that  have  been 
studied.  (See  also  2.14). 

Two  miles  downstream  of  the  narrows  that  mark  the  outlet 
of  Chilko  Lake,  the  river  is  very  wide  and  deep.  This  reach  of 


approximately  1  mile  can,  for  all  practical  purposes,  be  consid¬ 
ered  as  a  small  lake.  Between  the  two  lakes  the  Chilko  River  is 
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wide  and  slow  flowing.  There  are  a  few  islands  in  the  channel, 
but  one  could  hardly  call  the  river  braided.  This  section  of 

the  river  appears  to  be  stable.  No  traces  of  recent  river  eros¬ 

ion  could  be  seen.  A  measuring  reach  was  selected  between  the 
two  lakes.  It  will  be  described  in  2.14  (Fig*  41). 

After  leaving  the  lower  lake,  the  Chilko  River  enters  a 
valley  which  resembles  the  river  valleys  of  central  Alberta. 

It  is  between  600  and  2000  feet  wide,  and  100  to  200  feet  deep. 
Air  photos  show  terraces  and  abandoned  river  channels  at  various 
elevations  above  the  present  river  (Fig.  32).  A  comparison  of 

the  curvature  of  the  valley  sides  with  the  meander  pattern  of 

the  river,  suggests  that  a  larger  river  than  the  present  one 
has  done  the  major  part  of  the  valley-cutting  work.  Along  the 
first  ten  river  miles  only  a  few  places  of  active  river  eros¬ 
ion  could  be  detected  on  the  aerial  photographs.  Even  in  places 
where  the  river  impinges  on  a  steep  valley  side  there  appears 
to  be  little  erosive  activity.  Air  photos  show  trees  and  bushes 
growing  along  the  river  on  the  outside  of  such  curves. 

Approximately  5  miles  below  Chilko  Lake,  Lingfield  Creek 
joins  the  Chilko  River.  This  creek  is  almost  dry  during  the 
summer  months,  but  in  spring  it  might  have  a  discharge  of  200 
to  300  cfs.  It  is  steep  and  carries  a  heavy  load  of  coarse 
gravel  (picture  15). 

From  the  outlet  of  the  small  lake  (3  miles  below  Chilko 
Lake)  to  the  mouth  of  Lingfield  Creek,  the  Chilko  River  has 
one  distinct  channel.  Downstream  of  Lingfield  Creek  the  river 
is  partly  braided.  The  ratio  of  talweg  length  to  valley  length 
varies  from  1.2  to  1.4,  with  a  maximum  of  1.6  along  one  mile. 

Chilko  means  "clear  water”  in  the  language  of  the  Chilcotin 
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Indians.  (Chilcotin  -  people  of  the  clear  water.) 

The  Measuring  Reach.  The  lower  study  reach  on  the  Chilko 
River  begins  about  8  miles  below  Chilko  Lake  and  ends  3,600  ft. 
further  downstream  near  the  bridge  of  Henry's  Crossing  (Fig.  31 
and  Fig.  32).  The  reach  is  fairly  straight.  Along  most  of  the 
length  the  Chilko  has  only  one  channel,  but  there  is  a  distinct 
tendency  towards  braiding.  The  pictures  16  and  17  were  taken 
along  this  reach. 

Flow  Records.  The  Dominion  Water  Resources  gauging  stat¬ 
ion,  8MA2,  (Chilko  River  near  Tatla  Lake,  B.C.),is  located  at 
the  outlet  of  Chilko  Lake.  Thirty-four  years  of  records  exist 
for  this  station.  Fig.  3  shows  the  flood  frequency  curve  of 
this  station.  On  Fig.  33  the  mean  daily  discharges  of  the  last 
three  summers  have  been  plotted.  The  summer  of  1962  was  un¬ 
usually  wet  over  the  whole  Chilcotin  plateau.  Therefore  the 
July  discharges  of  the  above  gauging  station  have  been  increased 
by  5%  for  use  on  this  measuring  reach. 

Sections  and  Slope.  The  cross  sections  are  plotted  on 
Fig.  34  and  35.  Unfortunately  most  of  them  are  incomplete. 

Depth  was  estimated  with  equation  2.6.1.  The  general  shape  of 
the  bed  could  be  seen  clearly  from  the  banks.  The  section  at 
station  5  was  obtained  on  August  18,  1959,  by  engineers  of  the 
International  Pacific  Salmon  Fisheries  Commission  ( IPSFC) . 

They  spanned  the  river  with  a  cable  and  towed  an  echo  sounder 
across.  This  section  was  used  for  testing  equation  2.6.1. 

The  depth  calculated  from  this  equation  is  0.25  ft.  greater 
than  the  depth  obtained  by  dividing  the  area  of  the  section  by 
an  average  width,  (128  ft.).  At  the  bridge  the  river  has  been 
obstructed  by  material  that  was  pushed  into  it  from  the  left. 
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The  water  surface  slope  measurements  are  plotted  on  Fig. 

36.  On  Fig.  37  the  river  slope  (water  surface)  between  Cbilko 
Lake  and  Henry's  Crossing  is  shown.  It  has  been  measured  for 
tlie  B.C.  Hydro  and  Power  Authority  for  a  power  development 
study.  Water  surface  elevations  were  obtained  at  all  the  points 
for  which  the  talweg  length  is  shown  on  Fig.  37.  For  compari¬ 
son  the  slopes  that  were  measured  by  the  author  and  by  the  IPSFC 
are  also  shown. 

Table  7  gives  a  summary  of  the  channel  geometry. 

Bed  and  Bank  Material.  A  summary  of  the  bed  material 
samples  is  given  in  Table  2e  and  they  are  plotted  on  Fig.  38 
and  39.  On  Fig,  38  all  the  samples  that  were  collected  by  us¬ 
ing  some  grid  system  (walking  or  with  tape),  are  plotted  against 
frequency  of  exceedence  by  number.  On  Fig.  39  the  samples  that 
were  collected  over  a  certain  area  are  plotted  by  weight  and  by 
number.  Samples  2,  3,  and  5,  are  most  representative  of  the 
river  bed.  They  are  unfortunately  all  small  samples,  and  are 
therefore  unreliable. 

Fig.  40  shows  the  sieve  curves  of  two  gravel  samples  that 
were  taken  from  the  banks  of  this  reach. 

The  pictures  16  to  18,  ildusfrafe  t,he  1  drat  ion  and  the 
appearance,  of  several  grave  1-  samples . 

Regime  Considerations.  The  discharge  of  this  river  was 
unfortunately  quite  low  when  it  was  visited  first  in  July  1962. 
At  3,600  cfs  no  coarse  gravel  was  being  transported.  Had  trans¬ 
port  been  taking  place  it  would  undoubtedly  have  been  heard. 
Based  on  the  arguments  put  forward  in  2.7,  the  dominant  dis¬ 
charge  will  be  taken  as  5,600  cfs  which  corresponds  to  a  5  year 
flood.  The  dimensions  of  an  average  section  at  this  discharge 
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are  listed  in  Table  7.  Depth  has  been  calculated  from 
Manning's  equation. 

2 . 14  THE  CHILKO  RIVER  AT  THE  OUTLET  OF  CHILKO  LAKE 

A  general  description  of  the  Chilko  River  and  of  the 
flow  records  available  for  it  is  given  in  2.13. 

The  measuring  reach.  The  second  study  reach  on  the  Chilko 
River  lies  between  Chilko  Lake  and  the  small  lower  lake,  (see 
Fig.  31  and  Fig.  41).  This  reach  was  selected  because  the  bed 
pavement  appeared  to  be  finer  than  in  any  other  study  reach, 
and  because  the  International  Pacific  Salmon  Fisheries  Commis¬ 
sion  ( IPSFC)  had  surveyed  the  river  along  this  reach.  (The 
Chilko  Sockeyes  spawn  in  this  area),.  At  the  upstream  end  of 
the  measuring  reach  there  are  a  few  large  boulders  in  the  chan¬ 
nel.  At  low  stage  some  white  water  forms  there. 

Picturer,19  shows  this  reach. 

Sections  and  Slope.  The  cross  sections  are  shown  on  Fig. 
42.  The  main  channel  was  sectioned  by  echo  sounding.  The 
author  can  offer  no  explanation  for  the  wide  shallow  bench 
along  the  right  bank  of  the  river. 

The  slope  measurements  could  be  tied  into  the  slope  survey 
of  the  IPSFC.  All  the  water  surface  slope  measurements  are 
shown  on  Fig.  43.  The  IPSFC  slope  was  measured  with  stadia  at 
all  those  places  where  a  steep  river  bank  made  levelling  ted¬ 
ious.  This  should  explain  the  strange  "hump"  near  station  2, 
as  the  left  bank  is  steep  there. 

The  slope  and  section  measurements  are  summarized  in 
Table  8. 

Bed  Material.  The  results  of  the  bed  material  sampling 
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are  summarized  in  Table  2f  and  plotted  on  Fig.  44.  Samples  1 
and  2  were  taken  near  the  outside  edge  of  the  shallow  bench  on 
the  right  side  of  the  river.  Viewed  from  a  boat,  the  bed  mat¬ 
erial  of  the  main  channel  appeared  to  be  somewhat  coarser.  A 
4of  4  inches  will  be  assumed  for  that  part  of  the  channel. 
Practically  no  angular  stones  could  be  found,  but  in  each 
sample  about  one  quarter  of  the  stones  were  split,  probably  as 
a  result  of  frost  action. 

Regime  Considerations.  This  reach  lies  only  a  short  dis¬ 
tance  below  a  lake.  Therefore  a  dominant  discharge  of  low 
frequency  of  occurrence  seems  justified.  It  will  be  taken  as 
7,000  cfs,  which  corresponds  to  a  50  year  flood,  and  is  200  cfs 
greater  than  the  greatest  flood  on  record.  Using  Manning's 
equation,  the  dimensions  of  section  1  at  dominant  discharge 
have  been  calculated  #  The  results  are  listed  in  Table  8. 

2.15  THE  THOMPSON  RIVER  BELOW  KAMLOOPS  LAKE 

Engineers  of  the  B.C.  Water  Rights  Branch  have  recently 
been  investigating  the  possibilities  of  lowering  the  flood 
levels  in  Kamloops  City  by  improving  the  flow  conditions  at 
the  outlet  of  Kamloops  Lake.  A  large-scale  map  of  the  first 
4  miles  of  the  Thompson  was  prepared.  Fig.  45  is  a  copy  of 
the  key  plan  accompanying  this  map.  Seventeen  cross  sections 
were  obtained  by  echo  sounding.  The  water  surface  slope  was 
unfortunately  only  measured  at  low  stage  and  with  stadia. 

The  Measuring  Reach.  Along  the  first  mile  below  the  lake 
the  river  sections  show  extreme  variation  both  in  area  and 
depth  and  the  water  surface  slope  is  not  uniform.  Therefore 
only  the  reach  from  sections  7  to  17  (10,000  ft  along^of 

river), has  been  used  in  this  study.  The  channel  is  reasonably 


3  5 

straight  and  deeply  incised  (see  Fig.  46).  The  steepness  of 
the  banks  suggest  that  the  river  is  degrading  (see  picture  20). 
Bedrock  is  exposed  in  a  few  spots  on  the  left  side  of  the  chan¬ 
nel.  The  river  bed  seems  to  consist  of  coarse  gravel.  The 
possibility  of  bedrock  exposures  on  the  river  bed  cannot  be 
exc luded . 

Flow  Records.  The  Thompson  River  is  gauged  near  Spences 
Bridge,  some  60  miles  below  Kamloops  Lake.  (Thompson  River 
near  Spences  Bridge,  8LF51) .  Forty-seven  years  of  records  ex¬ 
ist  for  this  station.  In  1960  another  gauging  station  was  in¬ 
stalled  about  1  mile  below  Kamloops  Lake.  (Thompson  River  near 
Savona,  8LF33.)  Comparison  between  the  records  of  the  two 
stations  indicates  that  the  discharges  near  Savona  are  approx¬ 
imately  5%  smaller  than  near  Spences  Bridge. 

The  flood  frequency  curve  of  the  Thompson  River  near 
Spences  Bridge  is  shown  on  Fig.  3.  Based  on  this  curve  the 
fifty  year  flood  at  the  outlet  of  Kamloops  Lake  has  been  esti¬ 
mated  as  135,000  cfs  and  this  will  be  taken  as  the  dominant 
discharge.  The  flow  duration  curve  of  the  Thompson  River  near 
Spences  Bridge  is  shown  on  Fig.  47. 

Sections  and  Slope.  The  water  surface  slope  along  the 
measuring  reach  is  as  uniform  as  one  could  expect  from  stadia 
measurements  made  at  low  river  stage.  The  overall  slope  st.  7 
to  st .  17  should  be  reasonably  representative  of  the  river  slope 

at  high  stage.  It  is  0.72  parts  per  thousand. 

On  Fig.  48,  4  typical  sections  (9,  10,  15  and  16)  have  been 
plotted.  Section  15  is  similar  to  sections  8,  12  and  14. 

Section  9  is  similar  to  sections  10  and  11. 

Section  16  is  similar  to  section  13. 

Section  10  is  similar  to  sections  7,  11  and  17. 
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From  the  rating  curve  of  the  gauging  station  8LF33,  it  has 
been  estimated  that  the  water  surface  at  dominant  discharge  is 
17  feet  higher  than  the  water  surfaces  shown  on  Fig.  42.  On 
all  the  sections  the  vegetation  line  is  indicated.  The  water 
reaches  this  line  at  a  discharge  of  approximately  110,000  cfs, 
which  corresponds  to  the  five  year  flood.  A  summary  of  the 
cross  sections  at  dominant  discharge  is  shown  in  Table  9,  to¬ 
gether  with  some  other  relevant  information  about  this  measur¬ 
ing  reach. 

Bed  and  Bank  Material.  The  bed  pavement  was  sampled  on  a 
bar  near  section  10.  Equivalent  spherical  diameter  and  inter¬ 
mediate  axes  are  plotted  on  Fig.  49.  The  sample  is  also  listed 
in  Table  2.  The  author  is  unable  to  say  to  what  extent  the 
gravel  pavement  on  this  particular  bar  resembles  that  in  the 
deeper  parts  of  the  river.  The  B.C.  Water  Rights  Branch  party, 
who  surveyed  the  river,  believes  that  the  channel  is  paved  with 
this  type  of  material,  but  how  coarse  it  is  they  could  natur¬ 
ally  not  say.  Picture  20  shows  the  bar  on  which  the  gravel 
pavement  was  sampled  and  picture  21  is  a  close-up. 

Two  samples  were  taken  with  the  shovel;  one  of  gravel  from 
the  river  bank  near  station  4,l.b»  and  one  of  sand  from  a  pro¬ 
tected  little  beach  near  station  5,  r.b.  The  sieve  curves  are 
shown  on  Fig.  44. 
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CHAPTER  3  LABORATORY  WORK 
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3.2  EXPERIMENTAL  SETUP 

As  natural  a  gravel  as  possible  was  to  be  used  for  the 

& 

experiments,  but  the  equipment  available  in  the  hydraulics  lab¬ 
oratory  imposed  some  limiting  conditions  on  the  type  of  material 
that  could  be  employed.  From  the  experimental  data  of  Meyer- 
Peter  (unpublished)  it  was  concluded  that  in  order  to  obtain  a 
no-transport  equilibrium  at  sub  critical  flow,  and  at  a  dis¬ 
charge  of  1  to  2  cfs  (explained  below)  the  gravel  had  to  be 
seived  through  a  -|-inch  mesh.  The  laboratory  has  no  special 
dirty  water  system,  therefore, to  keep  the  main  water  system 
reasonably  clean,  the  fines  also  had  to  be  taken  out  of  the 
gravel . 

Samples  of  all  the  major  gravel  pits  in  the  vicinity  of 
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Edmonton  were  sieved  out.  Three  typical  sieve  curves  are  shown 
on  Fig.  51.  In  all  of  them  the  grain  sizes  between  0.7  mm  and 
5  mm  were  found  to  be  deficient.  This  was  not  desirable  for  the 
laboratory  work  but  it  had  to  be  accepted.  Ten  cubic  yards  of 
the  gravel  shown  in  Fig.  52  were  purchased.  In  the  laboratory 
this  material  was  then  sieved  through  a  -§--inch  mesh  and  washed 
on  a  50  mesh.  Sieve  curves  of  the  final  product  are  shown  in 
Fig.  53.  As  can  be  seen  the  washing  operation  did  not  meet  with 
success . 

A  glass  walled  flume,  60  feet  long,  4  feet  wide,  and  3  feet 
deep,  was  used  for  the  experiments.  It  was  set  at  a  slope  of 
0.5  percent  for  all  the  experiments.  In  the  first  test  equili¬ 
brium  was  to  be  approached  by  allowing  the  flow  to  build  up  a 
gravel  bed  in  the  flume.  For  this  run  a  short  inlet  chute  was 
installed  at  the  upstream  end  of  the  flume.  The  gravel  was  then 
fed  into  this  chute  by  means  of  a  conveyer  belt.  The  flume  was 
obstructed  at  intervals  of  a  few  feet  with  a  piece  of  plywood, 

2  inches  high,  to  start  the  building  up  process  on  the  smooth 
flume  bed.  The  flow  also  had  to  be  kept  away  from  the  smooth 
walls  in  order  to  simulate  natural  conditions.  This  was  done 
by  placing  pieces  of  plywood,  4  inches  wide,  along  the  sides. 
Some  coarse  gravel  was  also  put  into  the  flume  in  little  piles 
along  the  edges  before  the  first  run.  It  could  then  be  spread 
out  along  the  glass  walls  where  necessary. 

Discharge  was  measured  with  an  orifice  meter  in  the  over¬ 
head  pipe  between  the  pump  and  the  stilling  pond  at  the  head  of 
t he  f lume . 

As  the  general  aim  of  all  the  experiments  was  to  reproduce 
conditions  similar  to  those  encountered  in  the  field,  discharge 
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was  limited  to  one  or  possibly  two  cfs  by  the  width  of  the 
f 1 u me  . 

The  picture  22  shows  the  experimental  setup  before  the 
first  experiment. 

3.3  THE  EXPERIMENTS,  PROCEDURE  AND  RESULTS 
3.3.1  Test  1.  The  object  of  the  first  experiment  was  to 
approach  a  no-transport  equilibrium  state  by  building  up  a 
gravel  bed.  This  is  conparable  to  the  situation  on  the  river 
reach  ’’Quesnel  at  Lawless  Creek”,  where  the  river  bed  has  been 
built  up  by  the  Bullion  Mine  tailings. 

A  flow  of  1  cfs  was  started  in  the  flume  and  gravel  was 
fed  into  it  at  a  fairly  steady  rate  of  one  shovel  full  every 
five  seconds.  This  amounts  to  approximately  0.3  pound  per  sec¬ 
ond,  or  to  a  charge  of  0.5  percent.  Picture  23  shows  the  in¬ 
stallation  during  this  building-up  process.  Picture  24  was  tak¬ 
en  after  2  hours.  It  shows  how  the  building  up  proceeded  down 
the  flume  in  a  delta  like  manner,  with  a  steep  front,  2-3  inches 
high,  (visible  in  the  centre  of  the  picture).  Some  sand  that 
was  carried  in  suspension  past  this  delta  front  was  deposited 
further  down  the  flume  in  ripple  formation.  Flow  over  the 
gravel  bed  was  supercritical,  but  below  the  delta  front  it  was 
sub  critical. 

The  building  up  process  was  continued  for  5  hours>until 
the  delta  front  came  to  within  a  few  feet  of  the  tail-gate.  By 
then  the  upper  two  thirds  of  the  reach  were  mainly  gravel,  inter¬ 
spersed  with  sandy  patches  (picture  25).  The  sand  moved  along 
this  reach  as  sheet  flow,  practically  all  transport  taking  place 
close  to  the  bed.  Across  the  predominantly  gravel  areas  the 
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gravel  rolled  along  the  bottom,  coming  to  rest  occasionally. 

On  sandy  patches  the  gravel  grains  appeared  to  move  faster  and 
more  easily  and  rarely  came  to  rest.  The  deposits  in  the  upper 
parts  of  the  flume  were  firm.  They  had  a  gravel  skeleton  with 
sand  filling  the  interspaces. 

In  the  lower  third  of  the  flume  the  bed  was  predominantly 
sandy  (picture  26).  This  reach  was  somewhat  influenced  by  back¬ 
water  from  the  tailgate.  The  majority  of  the  gravel  grains  that 
entered  this  reach  rolled  across  it  and  were  deposited  at  the 
steep  delta  front.  The  deposits  along  this  reach  were  soft  and 
obviously  did  not  have  a  gravel  skeleton. 

Just  before  stopping  gravel  feeding,  the  water  surface 
slope  was  measured  at  five  foot  intervals,  in  the  centre  of  the 
flume  and  at  the  quarter  points.  Then  the  flume  was  immediately 
drained;  pictures  of  the  bed  were  taken  (25  and  26)  and  the  bed 
elevations  were  obtained  in  all  the  places  where  the  water  sur¬ 
face  elevations  had  been  measured.  Picture  26  shows  the  point- 
gauge  used  for  measuring  water  surface  and  bed  elevations. 

A  flow  of  1  cfs  was  then  resumed,  but  no  further  gravel 
was  supplied.  After  two  hours  of  degrading  transport  was  al¬ 
most  at  a  halt,  and  the  bed  was  distinctly  paved  with  the  coars¬ 
est  fractions  of  the  gravel.  The  flume  was  then  drained  and  bed 
elevations  were  obtained  as  before.  At  the  start  of  degrading 
the  tailgate  had  been  raised  0.13  ft.  to  provide  more  room  for 
the  deposition  of  material  at  the  lower  end  of  the  flume.  This 
shortened  the  significant  reach  (unaffected  by  entrance  and  exit 
conditions)  from  35  ft.  to  25  ft. 

To  ensure  that  a  state  of  equilibrium  had  been  observed,  a 
flow  of  1  cfs  was  again  started  and  left  running  for  twenty-four 
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hours.  At  that  time  absolutely  no  transport  could  be  seen. 

Water  surface  elevations  were  again  measured,  then  the  flume 
was  drained  and  bed  elevations  were  obtained.  No  change  great¬ 
er  than  0.02  ft.  had  taken  place  during  the  last  24  hours  of 
degrading.  Picture  27  shows  the  pavement  that  had  formed  on 
the  degraded  bed. 

The  author  expected  that  the  flow  would  cut  a  distinct 
channel  in  the  gravel  deposits  during  degradation,  but  this  did 
not  happen.  However  the  depth  became  somewhat  less  uniform 
across  the  channel.  The  variation  was  greatest  half  way  down 
the  flume  where  depth  was  0.18  ft.  on  the  right  side  and  0.10 
on  the  left  side.  Average  water  surface  and  bed  elevations  be¬ 
fore  and  after  degrading  are  shown  on  Fig.  54. 

3.3.2  Test  2.  As  the  first  test  had  not  produced  an  equili¬ 
brium  cross  section  resembling  a  natural  channel,  gravel  was 
lifted  into  the  flume  and  a  channel  was  constructed  on  the  de¬ 
graded  bed  of  Test  1  (picture  28).  The  channel  was  4  inches 
deep  and  15  inches  wide.  No  banks  were  built  in  the  lower  third 
of  the  flume.  A  flow  of  1  cfs  was  started  and  left  running  for 
24  hours.  After  that  time  no  transport  could  be  seen.  In  the 
upper  third  of  the  channel  a  gravel  river  with  gravel  banks  had 
formed  (picture  29).  Further  downstream  the  channel  had  wid¬ 
ened  and  the  depth  had  decreased,  until,  at  the  end  of  the  arti¬ 
ficial  banks,  the  flow  filled  again  the  whole  width  of  the  flume 
as  in  Test  1.  This  widening  was  caused  by  the  abrupt  ending  of 
the  artificially  installed  banks. 

Fig.  55  shows  the  average  water  surface  slope  after  equili¬ 
brium  has  become  established,  together  with  the  bed  elevations 
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along  the  centre  line  of  the  f luine  before  and  after  the  test. 

Two  cross  sections,  measured  15  and  20  f t  downstream  of  the 
flume  inlet,  are  also  shown. 

3.3.3  Test  3.  The  object  of  the  third  test  was  to  obtain  a 
no-transport  equilibrium  in  a  different  way  than  in  Test  1,  in 
order  to  establish^ whether  the  resulting  equilibrium  would  be 
identical  with  that  obtained  in  the  previous  test. 

The  gravel  was  shoveled  onto  a  very  steep  slope.  (Fig.  56). 

A  discharge  of  1  cfs  was  then  started.  This  experiment  is  com¬ 
parable  to  the  situation  than  exists  at  the  outlet  of  a  lake, 
which  is  dammed  by  glacial  debris. 

Initially  the  bed  degraded  very  rapidly  and  antidunes 
formed  in  the  lower  half  of  the  flume.  After  24  hours  all  trans¬ 
port  had  come  to  a  standstill.  Water  surface  and  bed  elevations 
were  then  measured  as  in  Test  1.  The  results  are  plotted  on 
Fig.  56. 

3.3.4  Test  4 .  A  final  test  was  made  by  running  a  discharge  of 
2.25  cfs  over  the  bed,  as  left  from  Test  3.  Again  no  transport 
could  be  seen  after  24  hours  of  flow,  when  the  flume  was  dewat¬ 
ered.  Water  surface  and  bed  elevations  were  measured  as  in 
Test  1  and  are  plotted  on  Fig.  56. 

During  this  test  the  baffle  piers  that  had  been  installed 
below  the  inlet  chute  were  washed  out.  Some  of  the  coarse  grav¬ 
el  that  had  been  used  to  pave  the  stilling  basin  was  moved  10 
to  20  ft  down  the  flume.  Therefore  the  steep  water  surface 
slope  at  the  upstream  end  of  the  flume  should  be  disregarded. 
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CHAPTER  4 _ WPLMAN'S  METHOD  OF  SAMPLING  GRAVELS 

AND  THE  ROUGHNESS  OF  CHANNELS 


During  the  course  of  the  field  work  the  method  of  gravel 
sampling  proposed  by  Wolman  (Ref.  11)  was  used  extensively. 

(See  paragraph  2.5  for  a  description  of  the  procedure.)  There 
are  two  different  ways  of  establishing  the  sampling  grid  and  the 
results  can  be  worked  out  in  several  fashions.  As  the  method  is 
relatively  new. no  particular  procedure  has  become  entrenched. 
Therefore  the  author  feels  justified  in  putting  forward  his 
opinion  about  the  proper  use  of  this  technique. 

4.1  THE  AIMS  OF  A  GRID  SAMPLING  TECHNIQUE 

It  is  accepted  standard  practice  to  characterize  the  grain 
size  distribution  of  non-cohesive  soils  by  taking  a  small  vol¬ 
ume,  sieving  it  out  in  square  mesh  sieves  and  plotting  the  re¬ 
sults  as  mesh  size  vs.  percent  of  total  dry  weight  passing 
through  that  mesh  (or  sometimes:  vs.  percent  of  total  dry 
weight  retained  by  that  mesh).  Since  this  kind  of  grain  size 
distribution  curve  has  been  used  for  practically  all  the  work 
done  on  rivers,  the  aim  of  Wolman's  grid  sampling  is  to  provide 
a  fast,  inexpensive  approximation  to  the  sieve  curve  of  the  top 
layer  of  a  gravel  deposit.  The  top  layer  is  usually  considered 
to  consist  of  all  the  material  between  the  surface  and  the  deep¬ 
est  exposed  particles.  (See  Fig.  57) f  (Ref.  11,  12,  14,  15). 

4.2  ESTABLISHING  A  GRID 

Wolman  suggested  two  methods  of  establishing  a  grid;  either 
to  use  actual  lines  and  to  collect  the  stones  vertically  below 
certain  points  (stretched  out  survey  tape,  checkerboard  arrange- 
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merit  of  wires,  etc.),  or  to  make  random  steps  picking  up  the 
stones  that  happen  to  fall  under  the  tip  of  one's  foot.  Qureshi 
only  used  the  second  method.  In  this  study  both  methods  have 
been  used.  Only  the  random  step  approach  is  feasible  for  samp¬ 
ling  under  water,  and  on  dry  ground  it  is  almost  impossible  to 
obtain  true  random  samples  by  stepping,  particularly  if  the  aver¬ 
age  diameter  of  the  stones  is  comparable  to  the  length  of  a  human 
f  oo  t . 

Unfortunately  the  two  different  grids  do  not  produce  simi¬ 
lar  results.  Table  2  and  Figs.  23  and  29  show  clearly  that  the 
grid  by  stepping  produces  samples  with  a  considerably  smaller 


). 


logarithmic  standard  deviation,  ( 2  6(,  ^ 


There  is  no  reason  to  doubt  that  the  samples  collected  by  both 
methods  are  reasonably  unbiased.  The  different  results  are 


caused  by  the  fact  that  the  operator's  foot  does  not  enter  into 


small  gaps  between  stones.  However,  if  a  grid  point  happens  to 
fall  above  a  gap, the.  usually  small,  stone  from  the  bottom  of 
the  gap  will  be  collected  (Fig.  58). 

The  curves  (log  b  or  log  vs.  percent  by  number  finer  than) 
produced  by  the  two  different  grids  do  not  differ  greatly  in  the 
coarser  part  of  a  sample.  As  this  range  of  the  curve  is  the 
most  important  one,  at  least  for  certain  applications  (see  4.5), 
it  seems  justified  to  continue  using  both  methods.  For  suffic¬ 
iently  large  samples  collected  with  a  fixed  grid  the  following 
statement  may  be  made:  "If  f  percent  (by  number)  of  the  stones 
in  a  sample  are  finer  than  ^  ,  then  f  percent  of  the  area  is 
covered  by  material  finer  than  4  This  statement  is  not  true 

for  samples  collected  by  pacing. 
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4.3  THE  INTERPRETATION  OF  SAMPLES 

4.3.1  The  Frequency.  Wolinan  takes  it  for  granted  that  the 
samples  have  to  be  plotted  against  frequency  by  number,  yet 
Qureshi  (Ref.  12)  used  frequency  by  weight  for  many  of  his 
curves.  The  following  points  will  show  that  the  use  of  frequ¬ 
ency  by  weight  is  very  time  consuming,  inaccurate,  subjective, 
and  also  does  not  produce  the  desired  results  as  stated  under 
4.1. 

1.  As  a  rough  approximation  to  a  natural  gravel,  one  can  assume 
a  material  consisting  of  only  two  grain  sizes,  D  =  1  inch  and 

D  =  y  inch,  each  size  accounting  for  50%  by  weight,  (Fig.  59). 

A  sample  of  N  stones  collected  from  this  material  by  using  a 
fixed  grid  will  probably  contain  roughly  ^  stones  of  each  grain 
size.  If  this  sample  is  then  analyzed  by  weight  the  median  dia¬ 
meter  and  even  the  lc-  will  be  1  inch,  while  a  plot  against  fre¬ 
quency  by  number  will  produce  a  curve  similar  to  the  sieve  curve 
of  the  material. 

2.  If  grain  size  is  plotted  against  frequency  by  weight  the 
resulting  curve  depends  almost  entirely  on  the  few  largest  stones 
in  the  sample.  On  river  beds  one  often  encounters  the  occasional 
boulder  that  does  not  really  belong  to  the  bed  material.  It 
might  be  an  erratic  boulder  that  fell  into  the  river  in  the 
courseof  meander  migration.  If  such  a  stone  happens  to  be  in  a 
sample  that  was  plotted  by  weight,  then  the  resulting  sieve 
curve  might  easily  indicate  grain  sizes  a  few  hundred  percent  in 
error.  To  overcome  this  difficulty  the  stones  greater  than  a 
certain  limit  would  have  to  be  rejected,  but  this  introduces  a 
subjective  factor  into  the  sampling  technique. 
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Ref.  12  suggests  a  rule  for  finding  the  maximum  size  that 
should  be  included  in  a  sample,  but  this  rule  is  so  complicated 
that  one  could  hardly  expect  a  field  party  to  apply  it  success¬ 
fully.  It  is  also  not  based  on  much  evidence  and  does  not  ex¬ 
clude  personal  judgement  completely. 

If  a  few  oversized  stones  happen  to  be  in  a  sample  that  is 
plotted  by  number  the  effect  is  almost  negligible.  Wolman’s 
paper  shows  how  well  samples  plotted  by  number  can  be  reproduced 
by  different  observers. 

3.  A  further  point  against  plotting  samples  by  weight  is  the 
amount  of  work  involved.  If  an  easily  measured  quantity  like 
the  intermediate  axis  is  plotted  against  frequency  by  number  the 
samples  can  be  analyzed  in  the  field  by  untrained  operators,  par¬ 
ticularly  if  a  simple  number  of  stones  like  49,  99  or  199  has 
been  collected.  (The  frequency  f  of  the  i-th  stone  in  a  sample 
of  n  stones  is  — -  .  The  use  of  n  +  1  instead  of  n  in  the  de¬ 
nominator  facilitates  plotting  on  probability  paper.) 

4.3.2  The  Linear  Grain  Size  Parameter.  The  size  of  the  inter¬ 
mediate  axis  determines  whether  a  grain  is  retained  or  passes 
through  a  given  round  hole  sieve.  The  same  is  almost  true  for 
square  mesh  sieves,  except  that  the  minor  axis  (c)  has  a  slight 
influence.  Wolman  (Ref. 11)  therefore  recommends  to  measure  only 
b  of  each  stone  and  to  use  it  as  the  linear  parameter  of  grain 
size.  Krumbein  and  Pettijohn  state  in  Ref.  16  that  the  equivalent 
spherical  diameter  and  b,  usually  lead  to  almost  identical 
results.  Qureshi  in  Ref.  12  used  for  some  curves  and  c  for 
others.  He  justifies  the  use  of  c  by  its  predominant  influence 
on  according  to  equation  2.5.2.  To  interpret  the  results  he 
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ol  tnined  by  plotting  c,  ho  uses  the  relation  <in.  =  1.55  c,.  .  Since 

'  r  f 

the  proportionality  constant  obviously  depends  on  the  type  of 
gravel  that  is  sampled,  and  since  c  by  itself  has  very  little 
in(  1  nonce  on  the  ordinary  volumetric  sieve  curve  of  a  gravel, 
the  author  does  not  consider  the  use  of  c  justified.  In  this 
study  all  the  major  gravel  samples  have  been  analyzed  by  using 
4-  and  b.  The  results  are  listed  in  Table  2  and  are  discussed 
separately  for  each  river  in  Chapter  2.  The  curves  are  re¬ 
produced  because  ot  the  work  involved  in  constructing  them.  The 
b  curves  look  very  similar.  A  typical  one  is  shown  on  Fig.  49. 
Table  2  indicates  that  the  logarithmic  standard  deviations  are 
almost  the  same  for  corresponding  b  and  cp  curves.  j  over 
b  ratios  at  constant  frequency  (  etc.  )  range  from  0.8 

b  &.<r 

to  1.03  and  stay  fairly  constant  for  a  given  river  reach.  Con¬ 
sidering  that  the  dp  ’s  of  equation  2.5.2  are  somewhat  small 
(Table  1),  a  good  average  ratio  is  0.95. 

Whether  p  or  b  should  be  given  preference  as  a  grain  size 
parameter  is  difficult  to  decide.  b  is  more  closely  related  to 
the  standard  sieve  curve.  From  the  point  of  view  of  a  worker 
interested  in  sediment  transport,  4*  offers  some  advantages  such 
as  its  close  relation  to  the  weight  of  a  stone  and  the  fact  that 
it  is  smaller  than  b  for  very  flat  stones.  Such  stones . rare ly 
ever  have  their  largest  face  exposed  to  the  full  force  of  the 
current.  Therefore  ^  might  characterize  the  exposed  surface 
more  adequately  than  b. 

For  future  work  the  author  recommends,  as  Wolman  did  prev¬ 
iously,  working  in  terms  of  b.  The  great  amount  of  labor  in¬ 
volved  in  calculating  cj}  is  not  justified  by  the  somewhat  doubt¬ 
ful  advantages  of  this  parameter.  All  three  axis  should  be 
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measured  on  a  few  typical  samples  to  have  some  record  of  the 
shape  of  the  material. 

4.4  COMPARISONS  BETWEEN  AREAL  AND  VOLUMETRIC  SAMPLES 

Wolman  concludes  on  the  basis  of  sieve  curves  (by  weight) 
for  volumetric  samples  and  frequency  (by  number)  plots  for  areal 
samples,  that  areal  grid  sampling  generally  produces  similar 
curves  but  indicates  larger  grain  sizes.  Qureshi  gives  sieve 
curves  and  areal  samples  that  were  taken  along  the  same  river 
reach.  The  sieve  curves  apply  to  the  top  layer  as  shown  in  Fig. 
57.  The  material  was  collected  above  water  at  low  river  stage. 
The  areal  samples  were  collected  from  below  water  by  stepping. 

The  average  based  on  approximately  15  sieve  curves  is  1.6 

inches,  D  r  is  1.9  inches.  The  average  cfV*  °f  the  eleven  samples 
for  which  or  c  curves  by  number  were  calculated  is  1.95 

inches  and  is  2.4  inches.  The  agreement  seems  very  close 

when  the  following  facts  are  taken  into  consideration:  The  vol¬ 
umetric  samples  were  taken  higher  up  in  the  channel,  where  the 
material  is  likely  to  be  finer,  and  the  samples  were  collected 
by  walking  rather  than  on  a  fixed  grid.  (See  4.2).  The  curves 
that  were  obtained  by  plotting  against  frequency  by  weight  give 
<p  values  several  times  larger  than  the  corresponding  D  values 
of  the  sieve  curves. 

4.5  THE  USE  OF  AREAL  SAMPLING  FOR  ESTIMATING  FRICTION  IN 

ROUGH  CHANNELS 

The  following  paragraph  is  based  mainly  on  Ref.  17  and 
Ref.  18. 

Friction  of  rough  turbulent  flow  is  commonly  described  by 
either  an  exponential  equation  of  the  form: 
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4.5.1 


or  by  a  logarithmic  equation  of  the  form: 
_  [ _  _  .  4?  ^ 

YF  ‘  ^  log  X 


4.5.2 


where 


R  =  hydraulic  radius 

-  constants  33R5 

f  =  friction  factor  =  - - 

\y vy 


k. 


-  Nikuradses  equivalent  sand  grain  roughness 


Since  it  is  well  established  that  the  energy  losses  in  fully 
turbulent  flow  over  rough  boundaries  are  closely  proportional 


to  vj  ,  equation  4.5.1  can  be  rewritten: 


IVl 


V. 


=  c. 


R  -C; 


k 


) 


4.5.1a 


v^  is  an  abbreviation  for  Wr?  and  is  usually  called  the 
friction  velocity. 

Equation  4.5.2  is  well  established  experimentally  for  pipes. 

The  coefficients  c.  and  c  t;.  depend  mainly  on  the  shape  of 
the  cross  sections.  Equation  4.5.1a  is  easier  to  use  than 


4.5.2  but  c  and  c  .. 


depend  on  the  shape  of  the  cross  sections 

R 

and  on  the  roughness  ratio  r~  . 

K  £ 


For  broad  channels  equation  4.5.2  was  given  by  Keulegan 


as  : 


VF 


1.5  +  2  log 


2  R 


4.5.2a 


This  can  also  be  written  as: 

V„ 


V, 


5.66  log  11.2 


11 

k. 


4.5.2b 


It  would  certainly  be  useful  if  ks  could  be  estimated  with 
Wolman's  sampling  method.  The  difficulty  lies  in  the  fact  that 
the  energy  losses  in  rivers  are  only  partly  the  result  of  grain 
roughness.  A  good  part  of  the  energy  loss  is  usually  due  to  bars, 
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bed  roughness  and  meandering. 

Of  the  seven  reaches  described  in  this  report  only  the 
reaches  "Quesnel  at  Lawless  Creek”  and  "Chilko  at  Henry's 
Crossing”  are  so  straight  and  uniform  that  a  correlation  be¬ 
tween  f  and  could  be  expected.  The  channel  of  the  Chilko 
River  contains  a  few  large  boulders  (see  picture  16),  but  none 
could  be  seen  in  the  Quesnel  River  (picture  2).  Values  of  ks 
have  been  calculated  from  equation  4.5.2a  for  all  the  river 
reaches  with  reasonably  uniform  flow  and  for  the  laboratory  test 
runs  1,  3  and  4,  (Table  11).  Instead  of  R,  d*  has  been  used 
t hr  oughout . 

Two  values  for  the  river  described  in  Ref.  14  are  also  in¬ 
cluded.  This  river,  the  Hasli-Aare,  is  a  gravel  river  that  has 
been  trained  between  rough  rigid  banks  some  80  years  ago.  On 
April  27,  1936>  it  was  at  low  stage  and  no  bed-load  transport 
took  place.  On  June  13,  1938,  it  was  in  full  flood  with  consid¬ 
erable  bed  load  transport. 

As  it  appeared  that  ,  b!/0  and  D?0  gave  the  best 
correlation  with  ks  (note  that  this  confirms  the  findings  of 
Meyer-Peter,  Ref.  19),  these  values  are  also  listed.  They  are 
based  on  areal  grid  samples  for  the  rivers  of  this  study  (Fig. 

10,  16,  29,  38  and  44  and  Table  2)  and  on  volumetric  surface 
samples  for  the  laboratory  data  (Fig.  53  and  Fig.  57)  and  for 
the  Hasli-Aare . 

The  agreement  between  k,  and  b?,,  or  D  /0  appears  very  good  for 
all  the  straight  uniform  channels.  This  again  indicates  that 
approximations  of  volumetric  sieve  curves  can  be  obtained  by 
areal  grid  sampling.  The  k  s  calculated  for  the  Cariboo  River 
at  Quesnel  Forks  is  somewhat  surprising.  The  badly  defined 
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slope  of  this  reach  might  be  responsible  (Fig.  15).  The  two 
k  5  values  calculated  for  the  Hasli-Aare  show  that  equation 
4.5.2a  does  not  hold  when  the  bed  is  in  motion.  For  the  reaches 
"Taseko  River  below  Taseko  Lake"  and  "Chilko  River  below  Chilko 
Lake"  the  fraction  of  the  total  energy  loss  that  is  spent  on 
grain  roughness  has  been  estimated  (F  ).  It  agrees  with  the 
values  suggested  in  Ref.  2. 

( !*• 

On  Fig.  6  0  771  has  been  plotted  against  for  the  follow- 

ing  data: 

(a)  The  straight  and  uniform  channels  of  Table  11.  !<<,  is  taken 

as  b  Cio  or  D^0  respectively. 

(b)  The  two  gravel  canals  described  in  Ref.  5.  ks  is  taken  as 
D  ?0  and  this  value  appears  to  be  based  on  a  volumetric 
surface  sample.  R  is  used  instead  of  d*  for  these  two 
points.  (Table  13). 

(c)  The  San  Luis  Valley  canals  described  by  Lane  in  Ref.  20. 

is  here  the  D  rlo  of  the  material  through  which  the 
canals  were  built.  It  is  possibly  somewhat  smaller  than 
the  D cj0  of  surface  samples  but  it  is  the  best  parameter 
available.  Lane  did  take  surface  samples,  but  not  as 
shown  in  Fig.  57.  (See  5.1  and  Table  13). 

The  curve  corresponding  to  equation  4.5.2b  and  two  straight 
lines  corresponding  to  equations  of  the  type  4.5.1a  are  also 
shown.  If  a  better  k5  value  had  been  available  for  Lane's 
canals,  there  would  probably  be  no  points  significantly  above 
the  line  of  equation  4.5.2b.  Some  scatter  below  this  line  has 
to  be  expected,  as  grain  roughness  of  the  channel  boundaries  is 
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ratios  of  less 


only  one  of  the  causes  of  energy  loss.  For 

'<s 

than  15  an  exponential  equation  of  the  4th-root  type  seems  to 
fit  the  data  best.  For  greater  ratios  a  line  corresponding  to 
Manning’s  equation  seems  to  fit  better.  (See  also  Ref.  17, 
p . 24  )  . 

4.6  AREAL  SAMPLING  WITHOUT  A  GRID 

A  few  samples  were  collected  by  picking  all  the  stones  with 
b  -  0.2  ft.  which  were  exposed  within  a  certain  small  area. 

This  corresponds  somewhat  to  the  ’’surface  samples”  of  Ref.  12,  21. 
Based  on  the  arguments  put  forward  in  4.3.1  such  samples  should 
be  plotted  against  frequency  by  weight,  to  give  something  com¬ 
parable  to  a  volumetric  sieve  curve.  Fig.  39  shows  plots  by 
weight  and  by  number.  It  is  interesting  to  note  that  the  fre¬ 
quency  curves  of  a  sample  taken  from  the  river  bed  pavement  and 
of  another  sample  taken  from  a  road  cut  through  the  river  bank 
plot  almost  identically.  This  sampling  technique  has  two  severe 
handicaps : 

1.  The  lower  limit  of  grain  sizes  to  be  considered  is  arbitrary, 
yet  without  such  a  limit  sampling  is  almost  impossible. 

2.  One  of  the  greatest  advantages  of  the  areal  sampling  with 
a  grid  is  lost,  namely  the  averaging  out  in  one  sample  of 
the  material  exposed  over  a  considerable  area. 
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CHAPTER  5  THE  EQUILIBRIUM  OF  GRAVEL  CHANNELS 


WITH  LOW  SEDIMENT  CHARGE 

5.1  OUTLINE  AND  ASSUMPFIONS 

In  this  chapter  an  attempt  will  be  made  to  find  the  rules  that 
correlate  the  stream  channel  geometry  of  a  special  type  of 
gravel  channel  with  discharge  and  sediment  size.  River  regime 
in  general  was  briefly  discussed  in  Chapter  1,  (See  particularly 
equation  1.3.2. 1).  Here  the  following  simplifying  assumptions 
will  be  made : 

1.  Straight  channels  with  uniform  flow. 

2.  Sediment  discharge  negligible  or  zero. 

3.  Discharge  constant  with  time. 

4.  Water  temperature  and  sediment  properties  other  than 
grain  size  distribution  unimportant. 

If  a  channel  meets  these  requirements  one  can  safely  assume 
that  the  dependent  variables  shown  on  the  right  side  of  equation 

1.3.2. 1  are  constant  with  time.  Of  these  variables  only  the 
slope  and  the  cross  section  as  defined  by  width  of  water  sur¬ 
face,  depth  and  area  will  be  of  interest  here. 

The  laboratory  data  meet  the  above  assumptions  almost  com¬ 
pletely,  but  the  seven  river  reaches  described  in  Chapter  2  only 
meet  them  partially.  Considerable  scatter  has  to  be  expected 
in  all  plots  because  of  this  and  because  of  incomplete  field 
data,  rough  estimates  of  dominant  discharge  and  uncertain  geo- 
morphological  conditions. 

The  main  conclusion  derived  from  Chapter  4,  was  that  grain 
sizes  obtained  from  grid  samples  and  worked  out  by  number  using 
4>  or  b  as  linear  parameters  are  very  similar  to  the  grain  sizes 
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obtained  by  sieving  volumetric  surface  samples  (whirl)  are 
analysed  by  weight)  of  the  same  areas.  It  was  also  concluded 
that  it  is  preferable  to  work  in  terms  of  the  coarser  diameters 
of  a  given  sample  (D^  ,  D  gs.  ,  ,  b6r  »  etc.)  because  in  this 

range  agreement  between  different  sampling  techniques  is  likely 
to  be  better.  These  two  conclusions  will  be  used  extensively 
in  this  chapter. 

The  U.S.B.R.  study  of  gravel  canals  (Lane,  Ref.  20)  will  be 
used  to  fill  the  wide  gap  in  discharge  and  grain  size  between 
the  laboratory  and  the  field  data.  These  canals  probably  car¬ 
ried  a  small  bed  material  load.  Lane  (Ref.  20  and  21)  states 
that  occasionally  coarse  material  entered  through  the  head  gates. 

In  the  hope  of  finding  some  characteristics  that  are  pecul¬ 
iar  to  gravel  rivers  with  low  sediment  discharge,  a  few  points 
for  channels  with  substantial  bed  load  will  be  included  in  all 
the  plots.  These  channels  are: 

1.  Laboratory  Test  1,  during  aggrading  (see  Table  10). 

2.  The  Aare  river  near  Brienzwiler  (Switzerland)  as  described 
in  Ref.  14.  A  discharge  of  6,720  cfs  (2  year  flood)  will 
be  taken  as  the  dominant  discharge  and  the  water  surface 
slope  and  average  cross  section  of  measuring  reach  1  (1,500 
feet  long)  will  be  accepted  as  representative.  At  a  dis¬ 
charge  of  5,000  cfs  a  bed  load  charge  of  90  ppm  was  meas¬ 
ured  . 

3.  Canals  12  and  13  of  the  field  data  of  Ref.  5.  They  are 
the  only  gravel  canals  in  this  set  of  observations.  They 
carry  a  heavy  sediment  charge  of  approximately  100  ppm, 
two  thirds  of  which  appear  to  be  bed  load. 
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The  hydraulic  data  of  (lie  rivers  and  canals  on  which  this 
chapter  is  based,  are  listed  in  Tables  12  and  13.  The  labora¬ 
tory  data  are  summarized  in  Table  10.  The  grain  size  distrib¬ 
ution  for  the  laboratory  tests  can  be  seen  on  Fig.  53.  Lane's 
bed  material  samples  are  comparable  to  the  areal  samples  without 
grid  (Table  2  and  Fig.  39),  which  the  author  collected  on  the 
reach  "Chilko  at  Henry's  Crossing".  To  make  his  results  compar¬ 
able  with  sieve  curves  of  the  top  layer  or  with  areal  grid  samples, 
the  author  adjusted  them  by  using  Figures  38  and  39. 

5.2  BRIEF  SURVEY  OF  REGIME  EQUATIONS 

The  concept  of  regime  originated  from  the  maintenance  of  un¬ 
lined  irrigation  canals  (Ref.  2,  Chapter  2).  Here  the  important 
independent  variables  are  essentially  constant  with  time,  the 
number  of  variables  is  reduced,  and  the  establishment  of  regime 
can  be  observed  without  great  difficulty.  From  the  observation 
of  irrigation  canals  several  competing  sets  of  equations  have 
been  developed  for  the  design  of  new  canals  (Ref.  2,  5-9).  The 
usual  approach  is  to  try  and  give  the  width,  depth  and  slope  of 
a  canal  in  terms  of  water  discharge.  The  great  majority  of  the 
regime  equations  make  W  proportional  to  Q‘  ,  d  or  R  proportional 

i  _  t  _  JL 

to Q  and  S  proportional  to  between  Q  ?  and  Q  ^  .  The  difference 
lies  in  the  way  in  which  the  proportionality  constants  are  cor¬ 
related  to  bed  load,  sediment  discharge  and  bed  and  bank  mater¬ 
ial.  The  disagreement  is  more  pronounced  in  the  form  of  the 
relationships  than  in  the  actual  result,  because  the  canal  data 
from  which  the  various  equations  were  derived  are  almost  ident¬ 
ical. 

The  U.  S.  Geological  Survey  has  made  an  extensive  investi- 


gat  ion  of  rivers  by  using  data  from  selected  gauging  sites 
(Ref.  10).  The  following  average  relations  were  derived  from 
measurements  made  at  different  gauging  stations  within  one 
drainage  basin.  Q  refers  to  discharges  of  equal  frequency  of 
exceedance . 
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This  seems  to  indicate  that  regime  equations  of  the  type  out¬ 
lined  above  can  be  applied  to  rivers. 

This  is  not  the  place  to  discuss  the  various  sets  of  regime 
equations.  However,  the  basic  relations  which  led  to  Blench's 
and  Lacey’s  equations  (Ref.  2,  Ref.  6-9)  have  to  be  outlined 
briefly,  because  their  applicability  to  the  gravel  channels  of 
this  study  will  be  investigated  in  paragraph  5.3. 

Lacey  found  that  vm  was  proportional  to  R2-  in  several 

Indian  canal  systems.  He  therefore  defined  a  silt  factor  f  as 

2 

f  =  0.73  ~  5.2.4 

The  0.73  was  included  to  make  f  =  1  for  a  well  known  set  of  canal 
data.  He  assumed  that  f  was  a  function  of  the  grain  size  of  the 
bed  material.  As  a  next  step  Lacey  concluded  that  A  f4-  ought  to 

Z  5 

be  a  function  of  vm  and  finds  A  f  =  3.8  vm  ,  which  leads  to 
vm  =  0.8  f3  Q1  5.2.5 

and 

P  =  2.67  Q  2  5.2 .6 

A 

where  P  stands  for  "wetted  perimeter”,  . 

The  plots  supporting  these  conclusions  are  shown  in  Ref.  9. 

In  Ref.  8  Lacey  deals  with  the  regime  slope.  On  the  basis 
of  Indian  canal  data  he  concludes  that 
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R'  S  =  constant  5.2.7 

for  a  given  type  of  bed  material  (in  Lacey’s  terms,  for  constant 
f  ).  This,  together  with  a  flow  formula  of  the  fourth  root  type 


(see  4.5.1a),  leads  to  his  regime  slope  formula 
v  =  16  (R2  S  )* 


5.2.8 


This  form  is  supposedly  independent  of  bed  and  side  material. 

2 

It  is  supported  by  a  plot  of  vm  vs  R  S,that  includes  a  wide 

range  of  bed  materials,  starting  with  gravel  rivers  of  the  Swiss 

Alps  and  ending  with  channels  on  silt. 

Blench  generalized  Lacey's  equations  to  differentiate  the 

effects  caused  by  the  sandy  beds  and  by  the  cohesive  sides.  He 

vJt 

calls  the  bed  factor.  It  is  supposed  to  depend  on  bed  mater¬ 

ial  and  bed  load  only.  A  side  factor  is  introduced? which  is 
meant  to  depend  on  the  erosion  resistance  of  the  sides  and  on 
the  tendency  of  suspended  load  to  deposit  on  the  sides.  A  cer- 
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tain  range  of  values  is  supposedly  possible  for  equilibri 


um 


VL 


conditions.  The  dynamical  explanation  for  given  in  Ref.  2 
applies  only  to  canals  with  smooth  cohesive  banks.  King's  re¬ 


gime  slope  formula,  quoted  by  Blench,  is 
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It  is  based  on  a  plot  of  — tt  vs.  — —  for  Indian  canals  with 

K  qdS  9 

steep  cohesive  banks  and  sand  beds  moving  in  dune  formation. 

For  constant  bed  and  side  factor  Blench's  regime  formulas 
satisfy  Lacey's  equation  5.2.7  in  terms  of  d  instead  of  R. 

When  Blench's  and  Lacey's  equations  are  written  as  design 
equations  to  give  R  or  d,  W  or  P  and  S  as  functions  of  Q,  silt 
factor,  bed  factor  and  side  factor  respectively,  the  exponent 
of  Q  is  found  to  be  the  same  in  corresponding  equations. 

Lane  (Ref.  20  and  21  approached  the  design  of  unlined 
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irrigation  canals  with  the  idea  that  a  canal  would  be  stable  as 
long  as  the  shear  force  exerted  by  the  flowing  water  on  the 
boundaries  did  not  exceed  a  certain  limiting  value.  This  criti¬ 
cal  or  limiting  shear  force  is  considered  to  depend  on  the 
diameter  of  the  grains  and  on  the  angle  at  which  the  material  is 
lying.  In  the  gravel  range  Lane’s  critical  tractive  force  curve 
is  well  supported  by  field  data. 

The  tractive  force  concept  cannot  be  reconciled  with  the 
regime  formulas  of  Lacey  or  Blench,  since  they  stipulate  that 

for  a  given  bed  and  side  material 
_1_ 

Rz  S  =  constant  5.2.7 

while  the  tractive  force  concept  leads  to 

R  S  =  constant  5.2.10 

for  a  broad  channel. 

5.3  EXAMINATION  OF  THE  DATA 

5.3.1  Regime  Formulas  of _ Lacey  and  Blench.  The  applicability 

of  the  regime  formulas,  outlined  in  paragraph  5.2,  to  a  certain 
set  of  regime  data  can  be  tested  in  various  ways.  Since  correl- 

2  _l 

ations  between  the  grain  size  of  the  bed  material  and  and  R"  S 
are  essential  for  Lacey's  and  Blench’s  formulas,  plots  of  grain 
size  versus  these  parameters  have  been  prepared  using  all  the 
channels  mentioned  under  5.1.  Grain  sizes  of  several  different 
frequencies  ("percent  finer  than")  were  used, but  none  of  the 

V  ~ 

plots  showed  any  functional  relationship.  The  use  of  instead 

of  did  not  produce  a  correlation  either.  Figures  61  and  62 

i  vtf 

show  the  plots  of  d** 2-  S  and  ~  against  the  "65  percent  finer 
than"  grain  sizes.  On  the  basis  of  such  plots  one  has  to  con¬ 
clude  that  regime  formulas  involving  silt  factor,  bed  factor  or 


side  factor  are  not  applicable  to  gravel  channels  of  the  type 
studied  here . 

Since  Lacey's  width  formula  (equation  5.2.6)  does  not  con¬ 
tain  the  silt  factor,  an  attempt  was  made  to  verify  it.  Fig. 

63  shows  a  plot  of  Q  vs  Ws  .  The  correlation  is  remarkable. 

The  best  fit  line  has  the  equation 

Ws  =  1.95  Q  *  5.3. 1.1 

Simple  explanations  can  be  put  forward  for  all  the  points  that 
fall  considerably  below  the  best  fit  line.  The  Hasli-Aare  was 
artificially  squeezed  between  rigid  banks  to  induce  it  to  lower 
its  bed.  The  Thompson  River  is  deeply  incised  and  a  few  rock 
outcrops  could  be  seen  along  its  sides.  Only  one  mile  below 
the  measuring  reach  it  widens  to  approximately  550  feet. 

Scatter  above  the  fitting  line  will  be  examined  in  5.3.3. 

5.3.2  Tractive  Force.  The  tractive  force  concept  was  tested  by 
plotting  grain  size  vs.  ^d*S.  Fairly  good  correlations  were  ob¬ 
tained  for  several  grain  size  frequencies  (percent  finer  than). 

On  Fig.  65  j'd^S  is  plotted  against  the  "50%  finer  than"  grain 
sizes.  The  fitting  line  of  this  plot  has  the  slope  of  one  upon 
one,  which  corresponds  to  Lane's  findings.  The  equation  of  this 
line  is 

D5o  =  0.38  y  d*S  5. 3.2.1 

The  best  correlation  was  obtained  by  plotting  the  "90%  finer 
than"  grain  sizes  against  ^d*S,  as  shown  on  Fig.  64.  The 
equation  of  this  fitting  line  is 

4 

D  =  0.7  (  ^d*S)T  5. 3. 2. 2 

The  scatter  on  both  plots  is  considerable  but  one  should  keep  in 
mind  that  it  is  impossible  to  characterize  a  natural  channel  pave- 
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ment  with  just  one  grain  size  parameter. 

The  following  interesting  conclusion  can  be  drawn  from  Fig, 
64  and  65: 

The  points  that  corresponds  to  channels  with  considerable 
bed  load  transport  plot  on  both  sides  of  the  fitting  line.  This 
throws  some  doubt  on  the  use  of  tractive  force  as  a  transport 
criterion.  (This  has  been  done  by  Meyer-Peter  in  Ref.  19.)  From 
results  derived  from  the  laboratory  work,  the  author  believes 
that  tractive  force  is  a  valid  erosion  criterion  for  paved, 
straight  and  uniform  gravel  channels.  If  d*S  increases  beyond 
a  critical  value  the  bed  starts  moving.  Once  the  bed  is  fully 
in  motion  the  bed  pavement  disappears  and  the  apparent  roughness 
of  the  channel  decreases.  With  it  the  tractive  force  may  de¬ 
crease  to  a  lower  value  than  the  original  critical  one.  The 
above  only  applies  to  natural  gravels  with  a  wide  range  of  grain 
sizes . 

5.3.3  The  Uniqueness  of  the  Equilibrium  of  Gravel  Channels. 

1.  A  few  points  on  Fig.  63  plot  considerably  above  the  fitting 
line  of  equation  5.  This  seems  to  indicate  that  the  w idth  of 
gravel  channels  in  low  transport  equilibrium  has  to  meet  only  a 
certain  minimum  condition  (given  by  equation  5. 3. 1.1).  Rivers 
with  considerable  suspended  load  build  banks  out  of  suspended 
material  and  keep  their  width  at  a  value  that  is  delicately  ad¬ 
justed  to  discharge  and  to  sediment  load. 

The  rivers  of  this  study  however,  have  very  little  suspended 
load  and  are  therefore  unable  to  build  banks.  Their  width  seems 
to  fit  equation  5.3.1. 1  closely,  whenever  the  present-day  river 
is  responsible  for  its  channel.  If  the  channel  contained  a 
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larger  river  at  some  earlier  time,  as  appears  to  be  the  case  at 
the  Chilko  Lake  outlet,  the  width  will  not  adjust  to  the  new 
discharge.  The  same  naturally  applies  to  gravel  canals  that 
have  been  built  wider  than  necessary.  As  a  design  equation 
5. 3. 1.1  should  probably  be  written  as 

Wy  -  1.95  5.3.3. 1 

2.  The  Chilko  River  below  Chilko  Lake  also  plots  considerably 
above  the  fitting  lines  of  Fig.  64  and  Fig.  65.  This  indicates 
that  equations  5.3.2. 1  and  2  should  have  the  equal  sign  replaced 
by  an  "equal  or  greater  than"  sign  (— ) ,  in  order  to  be  applicable 
to  all  gravel  channels.  The  Chilko  Lake  outlet  appears  to  be 
over-stable  (the  channel  could  take  a  considerably  larger  dis¬ 
charge  without  being  eroded)  and  not  in  regime. 

3.  The  main  aim  of  two  of  the  laboratory  tests  was  to  investi¬ 
gate  the  uniquenessj/of  slope  and  depth  for  a  given  d  is  char  ge  , 

bed  material  and  width.  The  equilibrium  was  approached  by  build¬ 
ing  up  a  bed  in  Test  lb  and  by  degrading  an  oversteepened  bed  in 
Test  3.  The  resulting  depth  and  slope  are  very  close  but  not 
identical.  In  channels  with  fixed  width,  the  way  in  which  equi¬ 
librium  is  reached  therefore  appears  to  have  only  a  minor  in¬ 
fluence  on  the  final  no-transport  equilibrium. 

5.4  REGIME  EQUATIONS  FOR  GRAVEL  CHANNELS  WITH  SMALL 

SEDIMENT  LOAD. 

With  the  simplifying  assumptions  made  for  this  chapter  (see 
5.1)  the  general  regime  function  (1.3. 2.1)  reduces  to: 

Ws  ,  d*,  S,  vw  =  f ns  (Q,  bed  material)  5.4.1 

or  to  : 

W,  d,  S,  vm  =  f ns  (Q,  bed  material)  5.4.1a 
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The  author  prefers  equation  5.4.1  because  Wg  is  easily  measured 
in  the  field.  It  will  be  assumed  that  the  bed  material  can  be 
characterized  with  one  linear  parameter  like  <^l0  ,  D  ,  etc. 

Four  independent  equations  are  needed  to  solve  5.4.1.  So 
far,  three  are  available;  namely  the  equation  of  continuity, 

Q  =  v  W  d*  5.4,2 

the  width  equation  5.3. 1.1,  and  the  tractive  force  equations 
5.3.2. 1  or  5. 3. 2. 2. 

When  the  flow  of  water  through  a  gravel  channel  carries 
practically  no  bed  load  it  should  be  possible  to  treat  it  as 
flow  in  rigid  boundaries  by  applying  the  friction  formulas  dis¬ 
cussed  in  4.5. 

Rivers  and  canals  flowing  through  sand  form  ripples  and 
dunes  on  their  beds.  The  channel  roughness  therefore  depends  on 
the  flow  conditions.  This  is  the  reason  why  the  flow  formulas 
of  Lacey  and  Blench  (5.2.8  and  5.2.9)  differ  from  the  rigid 
boundary  formulas  mentioned  in  4.5. 

In  the  field  and  in  the  laboratory  no  ripples  or  dunes 
could  be  detected.  Lane,  who  was  able  to  inspect  his  canals 
when  they  were  drained,  could  not  see  any  either  (Ref,  20). 

This  fact,  combined  with  the  findings  of  paragraph  4.5  (partic¬ 
ularly  Fig,  60),  support  the  conclusion  that  a  rigid  boundary 
flow  formula  with  the  grain  size  as  roughness  height  should  be 
used  as  a  regime  equation  for  gravel  channels  in  low  transport 
equilibrium.  An  exponential  flow  formula  will  be  used  here  be¬ 
cause  all  the  other  regime  equations  are  also  exponential.  The 
fourth  root  equation 


6.6  ( 


5.4,3 


v 


5 


appears  to  be  the  best  fit  to  the  points  of  Fig.  60.  To  test 
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the  assumptions  about  tractive  force  and  fourth  root  flow 
formula  one  can  combine  the  two  by  replacing  k  s  in  5.4.3  with 
Ds<7  or  D  90  of  equations  5.3.2. 1  or  5. 3. 2. 2  respectively. 

Using  5.3.2. 1  gives 


v„  ^  (d*2  S)^ 


and  5. 3. 2. 2  leads  to 


v„  -  (d*  '  S) 


5.4.4 


5.4.5 


These  two  equations  have  been  tested  on  Fig.  66  and  67.  Equation 
5.4.5  appears  to  be  better  supported  by  the  available  evidence. 
From  now  on  the  fitting  line 

vm  =  10(d*2S  S)6  5.4.5a 

will  be  used  instead  of  5. 3.2.2. 

The  points  that  do  not  fall  near  this  fitting  line  on  Fig. 

66  correspond  to  channels  with  considerable  bed  load  charge  or 
to  lake  outlets  which  are  probably  over  stable  (Chilko  and  Cari¬ 
boo).  Fig.  66  therefore  appears  to  be  a  useful  tool  for  testing 
the  no-transport  equilibrium  of  a  gravel  channel. 

The  four  equations  which  are  necessary  to  solve  the  simple 
regime  problem  of  5.4.1  are  now  available.  They  are: 

Q  =  v^  W5  d*  5.4.2 

Ws  =  1.95  5.3. 1.1 

vm  =  6 . 6  (  )'  5.4.3 

_1_ 

5.4.5a 


v. 


10(d*2’S  sr 


If  they  are  rewritten  to  give  the  dependent  variables  Ws  ,  d*, 

S  and  vfV1  in  terms  of  Q  and  ks  ,  the  following  equations  result: 
W  =  1.95  Q  '  5.3.1.  1 

C.j  -  O.  I 

d*  =  0.159  Q  ks  5.4.6 

*  O-d  O.  ZS 

S  =  0.1 172  Q  k^  5.4.7 


o.  i 


3 . 22  Q  ks 


O.  I 


v,( 


5.4 .8 


■ 

:0 


G4 


k5  is  the  Nikuradse  grain  size  roughness  measured  in  feet. 
According  to  the  conclusions  of  4.5  it  can  be  replaced  by  D <f0  , 

The  exponents  of  Q  are  equal  to  those  found  by  the  U.S.G.S. 

(Ref.  10  and  equations  5.2.1,  5.2.2,  and  5.2.3). 

If  Q  is  replaced  in  5.4.6  by  using  5.4.7  one  obtains  a  new 

fitting  line  for  Fig.  64.  This  equation  is 

± 

kff  =  0.818  (rd*S)3 

which  is  not  very  different  from  the  original  one. 

The  above  equations  do  not  apply  to  channels  with  fixed 
width  (flume  experiments,  trained  rivers,  overly  wide  canals). 
For  such  applications  5.4.6  and  5.4.7  have  to  be  replaced  by 

L \g  -0.1 

d*  =  0.268  Q  ks  Ws  5.4.9 

-0-2  o.£s  o,g 

S  =  0.071  Q  ks  W5  5.4 „ 10 

These  last  two  equations  could  be  checked  against  Meyer-Peter * s 
gravel  transport  experiments  (unpublished)  using  only  the  few 
ones  with  extremely  low  transport  rates.  The  agreement  is  rea¬ 
sonably  good  particularly  if  one  takes  into  consideration  that 
Meyer-Peter 1  s  experiments  were  made  by  adjusting  depth  and  dis¬ 
charge  to  obtain  a  certain  transport  rate  and  a  certain  water 
surface  slope.  The  sequence  of  cause  and  effect  in  a  river  is 
evidently  different. 

5.5  SOME  DIFFERENCES  BETWEEN  SAND  AND  GRAVEL  RIVERS 

The  regime  formulas  that  have  been  derived  in  5.4  differ 
appreciably  from  the  regime  formulas  of  Lacey  and  of  Blench. 

This  indicates  a  distinct  difference  between  the  mechanisms  of 
sand  and  gravel  transport,  as  there  can  be  little  doubt  that 
these  last  formulas  represent  the  regime  of  irrigation  canals 
with  sandy  beds  quite  well  (Ref.  5-9).  The  author's  time  and 
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data  are  insufficient  for  a  thorough  study  of  the  reasons  behind 
this  difference  but  he  would  like  to  point  out  a  few  pieces  of 
evidence  that  seem  to  support  his  belief  that  gravel  and  sand 
rivers  cannot  be  expected  to  obey  the  same  laws. 

5.5.1  The  Pavement  of  Gravel  Channels.  Sundborg  in  Ref.  22, 
points  out  how  the  relatively  rare  flume  experiments  using  mater¬ 
ial  with  a  wide  range  of  grain  sizes  have  come  to  different  con¬ 
clusions  depending  on  what  type  of  material  was  used.  Experiments 
with  materials  predominantly  in  the  sand  range  have  shown  that 
the  coarsest  grains,  if  not  coarser  than  approximately  6  mm, 
tend  to  move  first  and  most  easily  by  sliding  across  the  sandy 
bed.  Experiments  in  the  gravel  range  have  indicated  that  the 
finer  a  grain  is  the  more  likely  it  is  to  move. 

The  latter  implies  that  a  gravel  river  must  have  its  bed 
paved  with  the  coarsest  material  available, if  it  is  in  a  low 
transport  equilibrium.  An  easy  way  to  see  this  is  by  imagining 
a  gravel  river  with  a  bed  that  is  not  paved  and  that  consequently 
consists  of  a  wide  range  of  grain  sizes.  Such  a  river  could  only 
have  low  transport  rates  by  being  overstable  (unable  to  move  a 
major  portion  of  the  bed  material)  and  therefore  not  in  regime. 

If  it  were  in  regime  the  transport  rates  of  the  smaller  grain 
size  fractions  would  be  high.  The  fact  that  the  rivers  mentioned 
in  Chapter  2  and  the  laboratory  tests  showed  distinct  gravel 
pavements  is  an  indirect  proof  of  this  second  part  of  Sundborg's 
statement.  The  first  section  of  the  statement  is  supported  some¬ 
what  by  observations  made  during  the  aggrading  part  of  test  1. 
Sandy  patches  formed  on  the  bed.  Gravel  grains  moved  across  such 
patches  with  great  ease.  They  hardly  ever  came  to  rest  on  a  sandy 
patch. 
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A  second  implication  of  Sundborg's  statement  is  that  river 
reaches  with  bed  material  predominantly  in  the  2-4  mm  range 
should  not  exist.  This  material  would  move  with  greater  facility 
than  coarser  grains  through  the  gravel  reaches  of  a  river  and 
more  easily  than  sand  in  the  sandy  reaches.  Because  of  this  ease 
in  moving  along  the  river,  accumulation  of  2-4  mm  material  is  un¬ 
likely  to  occur.  Kuiper  (Ref.  23)  and  Halferdahl  (L.  B.  Halfer- 
dahl,  Geologist  with  the  Alberta  Research  Council,  oral  communi¬ 
cation)  have  investigated  several  of  the  larger  rivers  of  the 
Canadian  Prairies.  They  both  noticed  abrupt  changes  from  gravel 
to  sand. 

For  the  same  reasons  the  2-4  mm  grain  sizes  should  be 
deficient  in  river  deposits  in  the  areas  where  the  change  from 
gravel  to  sand  occurs.  The  gravel  samples  taken  in  the  Edmonton 
area  seem  to  confirm  this  (Fig.  51), 

5.5.2  Bed  Roughnesses.  The  second  important  difference  between 
sand  and  gravel  transport,  the  absence  of  ripples  and  dunes  in 
certain  gravel  channels,  has  already  been  mentioned  in  5.4.  Ref. 
15  describes  hydraulic  transport  experiments  with  gravel  of  a 
wide  grain  size  range.  At  high  transport  rates  dunes  formed, 
but  at  low  transport  rates  the  bed  stayed  flat  and  developed  a 

pavement.  Whether  dunes  form  seems  to  depend  on  the  roughness 

d* 

ratio  - —  and  on  the  transport  rate. 
k$ 

Two  of  the  main  formulas  of  paragraph  5.4  (tractive  force 
5.3.2. 1  and  the  fourth  root  flow  formula  5.4„3  with  grain  size  as 
roughness)  definitely  do  not  apply  to  channels  with  ripples  or 
dunes  on  the  bed.  It  is  therefore  unlikely  that  the  regime  of 
the  gravel  rivers  that  do  not  form  ripples  or  dunes  on  their  beds 
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can  be  described  by  formulas  having  the  same  form  as  equations 

5.4.6  to  5.4.10. 

5.6  TUB  DOWNSTREAM  EFFECTS  OF  A  DAM 

Rivers  usually  lose  their  sediment  load  in  the  ponds 
created  by  major  dams.  Downstream  of  such  a  dam  the  rivers  have 
to  adjust  to  the  enforced  low  sediment  transport  rates  and  pos¬ 
sibly  also  to  a  lower  dominant  discharge. 

Sand  rivers  adjust  by  changing  their  channel  cross  section 
and  their  slope  (Ref.  2,  10).  Gravel  rivers  have  the  added  pos¬ 
sibility  of  changing  the  characteristics  of  the  bed  material  by 
developing  a  pavement.  Instead  of  lowering  the  slope  such  a 
river  will  degrade,  by  picking  the  finer  grains  off  the  bed, 
until  the  bed  pavement  is  coarse  enough  to  resist  further  degrad¬ 
ing  (Ref.  15  is  a  study  of  the  formation  of  bed  pavements).  In 
terms  of  the  equations  developed  in  5.4  the  river  will  degrade 
until  ks  meets  the  requirements  of  equation  5.4.10  for  the 
existing  width  and  slope. 
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APPENDIX  A 


SYMBOLS  AND  ABBREVIATIONS 


Symbols 

a  large  ] 

b  intermediate  ■  axis  of  a  stone 
c  small  , 

9  equivalent  spherical  diameter 

The  indices  of  a,  b,  c  and  f  (b  ,  4b.)  mean, 
otherwise  specified,  ’’percent  finer  than”  bas 
grid  samples  worked  out  by  number. 

D  grain  size  based  on  sieve  curve  of  volumetric 
out 

worked/by  weight.  The  index  f  stands  for  ”pe 
than” . 


A 

d 

d* 
W  s 

w 

p 

R 

Q 

v 

Vm 

vs 

s 

V 

r 


area  of  channel  cross  section 

depth  in  the  central  parts  of  a  channel 

A 

water  surface  width 

A 

cf  " 

wetted  perimeter 
hydraulic  radius 
discharge 
velocity 

average  velocity  of  a  cross  section 
surface  velocity 
slope 

kinematic  viscosity 
specific  weight 


unless 
ed  on  areal 

sample 
rcent  finer 


of  water 


- 
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g  acceleration  of  gravity 
f  f  r  i  c  t  i  on  fact  or 
ky  Nikuradse  grain  size  roughness 
°F  degree  Fahrenheit 
V  volume 

All  dimensional  equations  are  in  the  foot-pound-second 
sys  tern  - 


Abbreviations 


ASCE 

American  Society  of  Civil  Engineers 

IAHR 

International  Association  for  Hydraul 

ICE 

Institution  of  Civil  Engineers  (G,B.) 

IPSFC 

International  Pacific  Salmon  Fisherie 

U.S.G. S, 

U.  S.  Geological  Survey 

U.S.B.R. 

U.  S.  Bureau  of  Reclamation 

Ref  . 

ref  erence ( s ) 

Fig  - 

f igure (s) 

Pr  oc  . 

proceedings 

Trans . 

transactions 

st  . 

station 

l.b. 

left  bank 

r  .  b . 

right  bank 

M .  C . 

main  channel 

sq .  m. 

square  mile 

Research 


Commission 


- 


U.uo  t  tnrrastn  I 
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APPENDIX  P, 

TABLES 

TABLE 

1 

COMPARISON  BETWEEN 

CALCULATED  AND 

MEASURED  WEIGHT 

OF  STONES 

River 

No .  of 
Stones 

Weight 
(grams ) 

Calculated 

(grams 

Cariboo  at  Q.F. 

16 

9231 

9560 

Chilko  at  Bridge 

8 

1227 

1180 

Taseko 

5 

10,440 

9270 

Taseko 

2 

14,700 

12,400 

Taseko 

9 

2537 

2310 

Taseko 

10 

3775 

3420 

Aare  (Ref. 14) 

1 

63,000 

43,500 

1 

40,000 

43,500 

1 

10,000 

6,500 

1 

15,000 

8,700 

1 

85,000 

78,000 

1 

17,000 

11,800 

1 

15,000 

17,400 

1 

9,000 

7,200 

1 

60,000 

67,400 

1 

17,000 

23,500 

1 

9,000 

9, 100 

1 

55,000 

82,600 

1 

12,000 

12,400 

1 

3,000 

1 , 700 

Aare , (Ref .  14)  , 
total  sample 

14 

410,000 

343,300 

Weight 
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SUMMARY  OF  GRAVEL  SAMPLES 
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TABLE  3 

QUESNEL  RIVER  AT  LAWLESS  CREEK 


-  Discharges  on  days  when  hydraulic  measurements  were  made: 

July  2  >  1962  14,800  cf s 

Sept. 20,  1962  5,340  cfs 

-  Water  temperature  at  high  stage:  47°  F 

-  Slope  measurements: 


Date_ Reac  h_ Slope 


Water 

July 

1-6 

6.53. 10"* 

surf  ace 

September 

5.68. 10_i 

Energy 

July 

1-5 

6.33. 10~B 

line 

September 

5. 16. 10~S 

-  Summary  of  cross 

section 

measurements 

(July  2 

,  62)  : 

Width  of 

Est  ir 

Cross 

Water 

ma  ted 

Section 

Measured 

St  n . 

Surf  ace 

Depth 

Area 

No. 

(ft) 

(ft) 

eft2) 

(ft  sH) 

(ft  s'1) 

1 

202 

8.1 

1340 

11 

2 

200 

8.5 

1360 

10.9 

12.3 

3 

216 

7.5 

1240 

11.9 

4 

150 

7.8 

1020 

11.6 

14.6  * 

5 

230 

5 . 6 

1190 

12 

*  M.C,  carrying  approx.  80%  of  total  discharge 


-  Average  section  at  dominant  discharge: 


Width  of  water  surface 

200 

f  t 

Dept  h 

8.0 

f  t 

Area 

1300 

f  t 

Dominant  discharge 

14,800 

cf  s 

Flood  Frequency: 

Recurrence  Interval 

Discharge 

(years) 

(cfs) 

2 

13,400 

5 

16,700 

10 

18,600 

20 

20,300 

■  '  r  . 
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TABLE  4 


TUB  CARIBOO 

)  RIVER 

AT  QUESNEL 

FORKS 

-  Discharges  on  days  when 

hydraulic  measurements  were  made: 

June  29,1962 

12,000 

cf  s 

July  1  ,1962 

10,500 

cf  s 

Sept.  17,1962 

2,810 

cf  s 

-  Water  temperature 

at  high  stage: 

46°  -  48°  F 

-  Slope  measurements: 

Date 

Reach 

Slope 

Energy  Line 

June  29 

1-4 

4.19. 10"5 

1-3 

4,62. 10~J 

Water  Surface 

June  29 

1-5 

3.71. 10’* 

1-3 

4.81. 10"* 

-  Summary  of  cross 

section 

measurements:  (June  29,  62) 

Stn.  Width  of 

Cross 

No.  Water 

Section 

Measured 

Surf  ace 

Depth 

Area 

V'rr,  V5 

(ft  ) 

(ft  ) 

(ft2  > 

(ft  s-' )  (ft  s">) 

1  165 

9 

1260 

9.5  11.5 

2  275 

6.5 

1150 

10.4 

3  192 

8 

1170 

10.3 

4  226 

8.5 

1544 

7.8  11.2 

-  Average  section  at  dominant  discharge: 

Width  of  water  surface  190 

Area  1250 

Depth  8.5 

Dominant  discharge  12,000 

f  t 
ft2 
f  t 
cf  s 

-  Flood  frequency: 

Recurrence  Interval  Discharge 

(years  (cfs) 

Depth  of 
Section 

Average 
(ft  ) 

2  12,700 

5  15,000 

10  17,000 

20  18,500 

8.5 

9.3 

' :  7 
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THE  CARIBOO  RIVER  AT  THE  OUTLET  OF  CARIBOO  LAKE 


-  Discharges  on  days  when  hydraulic  measurements  were  made: 


July 

8, 

1962 

6560 

cf  s 

July 

10, 

1962 

6710 

cf  s 

Sept . 

23  , 

1962 

2180 

cf  s 

-  Water  temperature  at 

high 

stage:  46° 

-  48°  F 

-  Slope  measurements: 

Date 

Reach 

Slope 

Energy  line 

July 

10  2-6 

1.64. 10_i 

Sept . 

23  2-6 

2.07. 10'5 

Water  surface 

July 

8  2-7 

1.92. 10" 1 

July 

10  2-7 

1.92 . 10~2 

Sept . 

23  2-7 

2.12.10 

-  Highest  discharge  in 

1962  : 

10,900  cf s 

at  a  gauge  reading 

2.5  ft.  higher  than  the  reading 

of  July  10. 

-  Summary  of  cross  section  measurements 

(July  10,  62)  : 

Width  of  Depth  of 

Water  Main 

Section 

Measured 

Stn.  Surface  Channel 

Area  d* 

v„,  vs 

No.  (ft)  (ft) 

(ft2)  (ft) 

(ft  s’1)  (ft  s  ) 

2  233  11.9 

1845  7.9 

3.65 

3  275  9.5 

1560  5.7 

4.3 

4  274  9 

1563  5.7 

4.3  7 

5  247  7 

1218  5.0 

5.5 

6  280  7.5 

1680  6.0 

4.0  7.7 

-  Average  section  at  dominant  discharge: 

Width  of  water  surface 

260 

f  t 

Area 

2800 

fta 

Depth 

13.0 

f  t 

Dominant  discharge 

17,000 

cf  s 

Bankfull  discharge 

8,000 

cf  s 

-  Flood  frequency: 

Recurrence  Interval 

Discha  r ge 

Estimated  Depth  of 

(ye  ar s ) 

(cfs) 

Average  Section  (ft) 

2 

11,300 

10.1 

5 

13,200 

11.1 

10 

15,000 

12 . 1 

20 

16,300 

12.8 

50 

17,600 

13.9 

100 

20,600 

14.9 

■  '  c  i;i:' ;)  hh-t 
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TABLE  6 


THE  TASEKO  RIVER  AT  THE  OUTLET  OF  TASEKO  LAKE 


Discharge  on  days  when  hydraulic  measurements  were  made 


Aug. 

13 

3500 

cf  s 

Aug. 

15 

3330 

cf  s 

Aug. 

17 

3330 

cf  s 

Aug. 

19 

3170 

cf  s 

Aug. 

21 

3550 

cf  s 

Aug. 

22 

3560 

cf  s 

Sept . 

8 

2500 

cf  s 

Water  temperature  at  high  stage:  approximately  48°  F 

Water  surface  slope: 

Date _ Reach _ Slope _ 

Aug.  13,  62  1-5  2.87. 1  Cf  * 

Aug.  21,  62  1-5  2.91.10  3 

Channel  dimensions  at  station  2,  Aug.  21,  62 


Width 

of 

water 

surf  ace 

135 

f  t 

Depth 

4.8 

f  t 

Area 

528 

ft2 

V  m 

6.7 

ft  s~' 

Measu 

red 

8 

ft  s" 

rage 

channel  d: 

Lmensions  at 

d  ominant 

discharge 

Wid  th 

of 

water 

surf  ace 

150 

f  t 

Area 

800 

ft  2 

Depth 

6.2 

f  t 

V  m 

8.0 

f  t 

Domin 

ant 

d ischarge 

6400 

cf  s 

Flood  frequency: 

Recurrence  Interval  Discharge 

( ye  a r  s ) (  cf  s  ) 


2 

5 

10 

20 


5400 

6100 

6500 

6900 
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TABLE  7 

CHTLKO  RIVER  AT  HENRY'S  CROSSING 


-  Discharge  on  days  when  hydraulic  measurements  were  made: 


Aug.  18 

1959 

3500 

cf  s 

July  16 

1962 

3420 

cf  s 

July  18 

1962 

3630 

cf  s 

Sept.  11 

1962 

2540 

cf  s  . 

- 

Water  temperature  at 

high  stage 

48°  -  50°  F 

- 

Water  surface 

slope : 

Date 

Reach 

Slope 

July  18, 

62 

0-5 

5 

.03. 10** 

Sept . 1 1 , 

62 

0-5 

5 

.08. 10'3 

- 

Summary  of  cross  section  measurements:  (July  18,  62) 

Width  of 

Esti- 

Cross 

Water 

mated 

Section 

Measured 

Stn . 

Surf  ace 

Depth 

Area 

v™  V  s 

Nc 

> . 

(ft) 

(ft) 

(ft1) 

(ft  s'1 )  (ft  s'1) 

0 

165 

4.1 

575 

6.3  9.6 

1 

135 

4.2 

460 

7.9 

2 

** 

147 

4.1 

525 

6.2 

3 

140 

4.7 

565 

6.4  7.3 

4 

155 

4‘.  2 

545 

6.8  8.5 

5 

* 

137 

4.7 

600 

5.9 

*  IPSFC  section 

**  Assuming  90%  of  discharge  is  in  main  channel 


-  Average  section  at  dominant  discharge: 


Width  of  water  surface 

150 

f  t 

Depth 

5.5 

f  t 

Area 

690 

ft* 

Dominant  discharge 

5600 

cf  s 

Flood  frequency: 

Recurrence  Interval 

Discha: 

( ye  a  r  s  ) 

(cf  s 

2 

4800 

5 

5600 

10 

6100 

20 

6500 

9?:  PI 
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CHILKO  RIVER  AT  THE  OUTLET  OF  CHILKO  LAKE 

-  Discharge  on  days  when  hydraulic  measurements  were  made: 

July  24,  1948  4600  cfs 

July  22,  1962  3340  cfs 

Sept  13,  1962  2370  cfs 


-  Water  temperature  at  high  stage: 

48°  -  50°  F 

-  Water  surface 

s lope : 

Date 

Reach 

Slope 

July  22 

,  62 

1-3 

0.92. 10~x 

Sept . 13 

,  62 

1-3 

0.94 . 10'3 

*  July  22 

,  48 

3700  ft 

1.19. 10'* 

*  IPSFC 

-  Summary  of  cross  secti 

on  measurements:  (July  22, 

1962) 

Width  of 

Depth  of 

Area  of 

Area  of 

Water 

Main 

Cross 

Main  Channel 

Main 

Stn.  Surface 

Channe 1 

Section 

Cross  Section 

Channe 1  * 

No .  (ft  ) 

(ft  ) 

(ft2> 

(ft*) 

(ft  s'1) 

1  272 

5.2 

1138 

988 

3.0 

2  390 

5.1 

1192 

735 

3.6 

3  263 

5.2 

1026 

812 

3.3 

*  assuming  80-90%  of  discharge  is  in  main  channel. 


-  Average  section  at  dominant  discharge: 


Width 

of 

water  surface 

270 

f  t 

Depth 

of 

main  channel 

7.5 

f  t 

Area 

of 

mai n  cha nne 1 

1413  * 

ft* 

Total 

area 

1760 

ft2 

Dominant 

dis  charge 

7000 

cf  s 

Bankf  ul 1 

di scharge 

5000 

cf  s 

*  carrying  80%  of  discharge. 


-  Flood  frequency: 

Recurrence  Interval  Discharge 

(years) (cf  s  ) 


2 

4800 

5 

5600 

10 

6100 

20 

6500 

50 

7000 

'■'3 
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TABLE  9 

THE  THOMPSON  RIVER  BELOW  KAMLOOPS  LAKE 


-  Discharge  on  days  when  hydraulic  measurements  were  made: 


Aug. 

17, 

1961 

Aug. 

18, 

1961 

Sept . 

28, 

1962 

23, 700  cf s 
24, 000  cfs 
18, 600  cf s 


-  Water  temperature  at  high  stage:  48°  -  50°  F 


-  Water  surface  slope: 

Aug . 17- 18 , 1961 ,  from  section  7  to  section  17:  0.72.10 

-  Summary  of  estimated  cross  section  dimensions  at  dominant 
discharge : 


Station 

Number 

Width  of 
Water  Surface 
(ft) 

Area 

(1000  ft2) 

Average 

Depth 

(ft) 

Maximum 

Depth 

(ft) 

V„ 

(ft  s-') 

7 

487 

12.4 

32.2 

39.3 

10.9 

8 

442 

10.1 

29.1 

31.7 

13.3 

9 

405 

11.0 

46.5 

49 

12.2 

10 

498 

11.6 

31.8 

42.2 

11.6 

11 

472 

12.6 

39 

43 

10.7 

12 

425 

10.5 

33.2 

37 

12.9 

13 

400 

11.0 

34.2 

44 

12.2 

14 

545 

14.4 

33.2 

38.5 

9.4 

15 

490 

12.3 

31 

33.5 

11 

16 

340 

11.2 

40 

56 

12 

17 

610 

18.2 

35.2 

45 

7.4 

-  Average  section  at  dominant 

Width  of  water  surface 

discharge : 

450 

ft 

Depth 

34 

f  t 

Area 

11,500 

ft* 

Dominant  discharge 

135, 000 

cf  s 

Bankfull  discharge 

110, 000 

cf  s 

-  Flood  frequency: 

Return  Period 

Discha 

rge 

(ye  ars) 

(1000 

cf  s  ) 

2  94 

5  110 

10  119 

20  126 

50  135 


" 
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TABLE  10 

SUMMARY  OF  LABORATORY  EXPERIMENTS 


Test 

Dis¬ 
charge 
(cf  s) 

Average 
Water  Sur¬ 
face  Slope 
(%) 

Average 

Depth 

(ft) 

Width 

(ft) 

(ft  s") 

Measured 

V. 

(ft  s'  ) 

agr . 

1 

1.31 

0.10 

3  o  4 

2.80 

2.92 

1 

degr . 

1 

1.20 

0.14 

3.4 

2.18 

2.85 

2 

1 

0.99 

0.22 

2.6 

2. 10 

3.33 

3 

1 

1.24 

0. 15 

3.0 

2.27 

4 

2.25 

0 . 625 

0.26 

3  o  4 

2.57 

3  o  34 
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TABLE  11 

EQUIVALENT  SAND  GRAIN  ROUGHNESS 

k  5.  <Pic  b  q0  Dc(0 


Channe 1 

Date 

(ft) 

(ft) 

(f  t ) 

(ft) 

F 

Quesnel  at 

Lawless  Creek 

July  2 , 62 
Sept . 20 , 62 

1.24 

1 0  35 

1.08 

1.25 

Chilko  at 

Henry's  Crossing 

July  18,62 
Sept  oil,  62 

1 .  15 
1.17 

0.83 

0.83 

Laboratory 

Test  1 

Test  2 

Test  3 

0.031 

0.037 

0.030 

0.033 

Hasli  Aare  at 

Br ienzwiler 

Apr .27,36 
June  13,38 

0.61 

1 0  50 

0.59 

Cariboo  at 

Quesnel  forks 

June  29,62 

1.15 

1.4 

1.55 

Taseko 

Aug .21,62 

3  0  18 

0.95 

1.00 

0.47 

Chilko  at 

Chilko  Lake 

July  22 , 62 

1.55 

0.55 

0.59 

0.62 

-k 


is  calculated  from  equation  4,502a,  using  d*  instead  of  R 
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TABLE  13 


HYDRAULIC  DATA  OF  GRAVEL  CANALS  AT  FULL  SUPPLY 


(1) 

(3) 

(3) 

Q 

Ws 

d 

d  * 

A 

S  „ 

Dm 

Dfo.v 

D  „ 

Cana  1 

cf  s 

ft 

f  t 

ft 

ft1 2 3 4 

ft  s 

.  10 

in 

in 

in 

Lane 

(Ref. 

20) 

#1 

1500 

73 

4.87 

3.50 

255 

5.88 

2.80 

5.0 

2.69 

2.32 

tl 

#4 

768 

48 

3.11 

2.46 

118 

6.53 

3.59 

2.55 

2.49 

2.15 

It 

#5 

448 

40 

2.5 

1.95 

77 

5.82 

3.68 

2.65 

1.78 

1.51 

n 

#6 

159 

21.7 

1.88 

1.60 

34.6 

4.59 

2.95 

2.20 

1.38 

1.20 

M 

#7 

95.6 

15.9 

1.73 

1.38 

21.9 

4.36 

2.90 

1.30 

1.35 

1.15 

It 

#8 

46 

19.2 

0.96 

0.80 

15.3 

3.00 

3.16 

1.45 

1.28 

1.10 

it 

#10 

16.6 

11.1 

0.60 

0.515 

5.7 

2.90 

9.65 

2.50 

2.1 

1.80 

tt 

#11 

203 

32.3 

1.88 

1.62 

52.4 

3.88 

2.35 

1.70 

1.58 

1.36 

tl 

#12 

128 

21.9 

1.77 

1.46 

32.0 

4.00 

2.43 

1.45 

1 . 12 

0.96 

tt 

#14 

110 

21.4 

2.00 

1.56 

33.4 

3.29 

1.39 

0 . 66 

0.66 

0.56 

tt 

#15 

477 

39.4 

3.05 

2.50 

99 

4.84 

1.99 

2.45 

1.64 

1.42 

u 

#17 

531 

41 

2.60 

2.35 

96 

5.51 

2.74 

1.90 

1.25 

1.08 

Simons  (Ref. 5) 

(2) 

D  g.~ 

(4) 

(2) 

#12 

883 

60 

7.88 

6.52 

391 

2.26 

0.181 

0.98 

0.63 

0.276 

#13 

751 

60 

5.73 

4.87 

292 

2.57 

0.166 

0.79 

0.55 

0.300 

(1)  Based  on  sieve  analysis  of  material  through  which  the 
canals  were  built. 

(2)  Based  on  sieve  analysis  of  bed  material. 

(3)  Estimated  from  Lane's  surface  samples. 

(4)  Estimated  from  D  and  . 
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iUESNEL  HIVEH  AT  LAWLESS  CREEK 


Picture  1  Looking  down  from  road  between  Quesnel  Forks 
and  Likely  on  to  the  upper  part  of  the  meas¬ 
uring  reach. 

July  3,  1962. 


Picture  2  As  picture  1  but  on  Sept.  20,  1962,  at  low 
river  stage. 
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Picture  3  Area  of  gravel  sample  l}.,  July  4,  1962. 


cariboo  river  at  ^uesnel  porks 


Picture  4  Looking  over  the  measuring  reach  from  st.4, 
l.b.,  on  June  30,  1962. 
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Picture  5  Looking  upstream  from  st.3  on  Sept.  l8,  1962. 


Picture  b  Tape  stretched  out  for  gravel  sample  1,  Sept. 
19,  1962. 


Picture  7 


Road  cut  near  s t * I4. *  The  material  for  sieve 
curve  1  was  taken  here. 


t 


GARIBOC  RIVER  AT  OUTLET  OF  CARIBOO  LAKE 


Picture  0  Looking  downstream  from  st.5  on  July  10,  1962 


1  57 


Picture  9  Typical  bank  at  high  stage. 


Picture  1C  Place  of  gravel  sample  5  near  st.3,  on  bar 
in  centre  of  river. 


Picture  11  Gravel  sample  3,  at  waterline,  st.3,  r.b. 


TASEKO  RIVER  BELOW  TASEKO  LaKE 


Picture  12  Looking  downstream  along  the  measuring 

reach  from  a  place  between  st.l  and  st.  2, 
r.b.  Gravel  sample  5  was  taken  here.  Aug. 
19,  1962. 


Picture  13-  Looking  downstream,  from  st.  I4.,  l„b.  Gravel 
sample  3  was  taken  on  bar  in  foreground, 
Aug,  21,  1962. 


Picture  14  Obstructed  part  of  the  river  at  the  upper 
end  of  the  measuring  reach.  Tape  in  place 
for  gravel  sample  2.  Aug.  19,  19o2. 


CHILKO  RIVER  AT  HENRY 


CROSSING 


Picture  1$  Lingfield  Creek  in  September. 


Picture  lb  Looking  upstream  from  near  st.  $,  l.b. 

Gravel  sample  5  was  taken  from  the  almost 
exposed  bar  in  the  river.  Sept,  llj.*  19o2. 
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Picture  17  Location  of  gravel  sample  6  at  st.  1+.,  r.b. 
Sept.  12,  1962, 


Picture  l8  Material  exposed  in  roadcut  through  river 
bank.  (Sieve  curve  2,  gravel  sample  7) 


CMILKO  RIVER  AT  THE  OUTLET  OF  CHILKO  L 
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Picture  19  Looking  downstream  from  at.  2,  r.b.  on  July 
23,  1962. 


THOMPSON  RIVER  AT  THE  OUTLET  OF  KAMLOOPS  LAKE 


Picture  20  Looking  upstream  at  the  gravel  bar  on  which 
the  gravel  sample  of  Fig.  43  and  of  picture 
21  were  taken. 


Picture  21  Bed  pavement  on  bar  near  st.  10,  r.b. 


LABORATORY  WORK 


Picture  22  The  flume  before  the  fir3t  test 


104 


Picture  2k  The  dewatered  flume  after  2  h  of  aggrading 
(test  1) 


Picture  23  The  upper  part  of  the  flume  during  test  1, 
aggrading . 
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Picture  25 


The  bed  near  the  upstream  end  of  the  flume 
after  aggrading,  (test  1) 


Picture  26  Dewatered  flume  after  aggrading .( test  1). 

Point  gauge  set  for  measuring  bed  elevation. 


Picture  27  Close  up  of  bed  pavement.  Plow  was  left  to 
right . 


Picture  28  Artificial  channel  of  test  2. 
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